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Executive Summary 


This report summarizes research completed under Contract NAS 1-1 8239 to the 
University of Chicago. The effort has involved application of photochemical models to 
interpret measurements made by the Halogen Occultation Experiment (HALOE) and the use 
of HALOE data to examine stratospheric chemical processes which influence the Earth’s 
ozone layer. This work constitutes the Ph.D. research of Loretta J. Mickley at the 
University of Chicago. 

This research has emphasized the partitioning of members of the odd nitrogen (NO y ) 
and odd chlorine (Cl y ) families in the stratosphere. Use of two tracer molecules measured 
by HALOE, CH 4 and HF, enabled us to follow the partitioning of the families over time 
scales ranging from weeks to years. In one study (Project 1) we determined the effect of 
aerosol from the Mount Pinatubo volcano on NO x (NO and NO z ) and ozone at altitudes 
which are not assessable by research aircraft. In a second study (Project 2), we 
investigated the effect of the NO and N0 2 on Cl y partitioning during Antarctic spring, after 
significant ozone depletion had taken place. In addition to these projects, we also worked 
with scientists at NASA/Langley to validate HALOE measurements of NO and N0 2 . For 
this work, we compared the instrument’s measurements with observations from other 
sensors and with model-based predictions. A summary of the two major projects appears 
below. 

Project 1 : Response of Stratospheric Odd Nitrogen and Ozone at 17 mb to Mount Pinatubo 
Aerosol in the Southern Mid-Latitudes 

Observations of stratospheric NO, N0 2 , and 0 3 from HALOE were examined over the 
Southern mid-latitudes (35-45 degrees South) in the years following the eruption of Mount 
Pinatubo, from early 1992 through mid- 1995. The tracers HF and CH 4 , also observed by 
HALOE, were used to distinguish between the effects of transport on the distribution of 
ozone and total reactive nitrogen (NO y ) and the chemical effects of aerosol on NO y 



partitioning and ozone loss processes. Mixing ratios of HF were adjusted to take into 
account its increasing burden due to emission of chlorofluorocarbons at the Earth’s surface. 

Results show that at 17 mb (about 28 km), in parcels containing similar values of HF, 
the abundance of NO+N0 2 (NO,) almost doubled during the first two years of the 
observation period. Since model results show that HF correlates with NO y in this part of 
the stratosphere, the increase must be due to changes in the partitioning of NO y . Data on 
aerosol surface area from SAGE II show a sharp decline during 1992 at this altitude, 
followed by a leveling off. Thus the increase in NO, may be linked to a decline in the rate 
of hydrolysis of N 2 O s on sulfate aerosol as the Pinatubo aerosol gradually settled out of the 
stratosphere. Similar results were obtained when CH 4 was used as a tracer. 

Again using HF and CH 4 as tracers, results show that ozone abundances at 17 mb 

declined by about 10 percent during the first year of the observation period. The trend in 

ozone is opposite in sign to that observed by instruments at lower altitudes (18-22 km) 

following the Mount Pinatubo eruption. At the lower altitudes, NO, acts as a brake on the 

active chlorine and hydrogen cycles that destroy odd oxygen, and a return to pre-Pinatubo 

values of NO, resulted in an increase in ozone amounts. The ozone trend observed at 

higher altitudes is consistent with model calculations which show that above about 26 km, 

where NO, mixing ratios begin increasing rapidly with height, the NO, catalytic cycle 

dominates all other ozone loss processes. Thus the data suggest that in 1991 and 1992 at 

17 mb, the loss in NO, due to heterogeneous chemistry on Pinatubo aerosols may have 

«• 

temporarily elevated levels of ozone. 

Project 2 : Evolution of Chlorine and Nitrogen Species in the Lower Stratosphere during 
Antarctic Spring 

Observations of 0 3 , HCl, NO, and NO z from HALOE provided a means to investigate 
chemical change in the lower stratosphere over Antarctica during 23 days in October of 
1992. Two long-lived species also observed by HALOE, HF and CH 4 , were used as 



tracers to identify a series of air parcels having similar Cl y and NO y abundances and similar 
temperature histories. Most of the parcels chosen using tracer analysis show low ozone 
mixing ratios, less than 1 ppm on the 480 Kelvin surface (about 20 km). The low 0 3 
content of these parcels suggests that ozone depletion had already taken place, during the 
month of September. HC1 mixing ratios in these parcels more than doubled from early to 
late October, and the NO and N0 2 abundances rose sharply, by a factor of five for NO and 
a factor of three to four for N0 2 . 

The trends observed in NO, N0 2 , and HC1 agree with model calculations which show 
that the formation of HC1 proceeds quickly when ozone levels fall so low that (1) the rate of 
the reaction Cl + 0 3 — > CIO + 0 2 slows and (2) the rate of the reaction NO + CIO — > N0 2 
+ Cl becomes faster than the rate of the competing reaction NO + 0 3 --> N0 2 + 0 2 . Under 
these conditions, Cl increases at the expense of CIO, and HC1 is formed via the reaction Cl 
+ CH 4 — > HC1 + CH 3 . Stratospheric chlorine is thus shifted from reactive species to the 
long-lived reservoir molecule, HC1. The repartitioning of the active chlorine family in 
favor of HC1 may limit the extent of ozone depletion during spring, and may have 
importance as the polar vortex breaks up and parcels are transported toward mid-latitudes. 
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CHAPTER 1 


INTRODUCTION 

Ozone in the stratosphere is of major importance to life on earth. The gas 
absorbs strongly in the ultraviolet region of the spectrum (200 to 320 nm), light 
that in large enough doses can severely damage plant and animal life by breaking 
chemical bonds. By filtering out these wavelengths of sunlight, the stratospheric 
ozone layer protects the biosphere. The ozone layer also plays a role in the radiative 
balance of the atmosphere. Absorption of incoming sunlight by ozone prevents a 
fraction of that light from reaching the earth’s surface. The absorption also warms 
the stratosphere, influencing temperatures at the tropopause. 

In recent years, disturbing trends in stratospheric ozone have been observed 
over the mid-latitudes and polar regions. The first sign of a perturbation to the 
ozone layer was detected during spring over Antarctica (Farman et al., 1985). There 
ground-based instruments measuring total column amounts showed large losses of 
total ozone each October beginning in the mid-1970s. The Antarctic ozone hole, 
monitored now by satellites, balloons and aircraft as well as ground-based instru- 
ments, has occurred each spring since then, with especially severe depletions taking 
place in spring 1992 and 1993 (Hofmann et al., 1989a; Waters et al., 1993; WMO, 
1995). The data show losses in total ozone of 50-60 per cent over a region of 20 
million square kilometers in area, or from the pole to about 65 S. Nearly complete 
loss of ozone has been observed in the 14-17 km layer of the stratosphere (Hofmann 
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et al., 1994a). Similar, though not as severe, losses of ozone have been observed over 
the Artie in spring (Profitt et al., 1990; Browell et al., 1990; Waters et al., 1993). 

Beginning in the mid-1980s, long-term trends in ozone loss over the mid- 
latitudes have also been reported. Data from satellites, ozonesondes, and ground- 
based instruments all show trends in ozone loss of a few per cent per decade, with 
greater loss occurring in the Northern Hemisphere (Reinsel et al., 1984; Stolarski 
et al., 1992; Reinsel et al., 1994; WMO, 1995). These trends in ozone appear to 
have begun in the early to mid-1970s. Since total column ozone has been measured 
by ground-based spectrometers since the 1920s, the loss trends can be determined 
with some confidence. Briefly, the data show recent losses of about 4 to 5 per 
cent per decade over the northern mid-latitudes, and 3 to 4 per cent per decade 
over the southern mid-latitudes. In the Northern Hemisphere, losses are greatest 
in late winter through early spring. Data from the Stratospheric Aerosol and Gas 
Experiment (SAGE) and from a network of ozonesondes show that most of the loss 
over mid-latitudes occurs in the 15-20 km region (McCormick et al., 1992; Logan, 
1994). Over the tropics, ozone appears to have declined only a fraction of a per cent 
per decade (WMO, 1995). 

Recent work has been done to try to detect trends in ultraviolet flux at the 
earth’s surface that are consistent with these trends in ozone. Over mid-latitudes, 
problems in instrument calibration together with the challenge of removing the ef- 
fects of clouds, aerosol, and pollution have made it difficult to determine trends 
in ultraviolet flux due solely to trends in stratospheric ozone. (See, for example, 
Frederick et al., 1992.) Lubin and Jensen (1995) show that up till now the large 
natural variability in ultraviolet flux due to variable cloudiness probably swamps 



3 


any change in the irradiance due to loss of stratospheric ozone. Frederick and Erlick 
(1995) have reached a similar conclusion. However, changes in summertime ultravi- 
olet radiation due to changes in the ozone layer may in fact become significant by 
the end of this century (Lubin and Jensen, 1995). Over the southern high latitudes 
in spring, enhanced flux of ultraviolet sunlight has already been reported (Frederick 
and Alberts, 1991). 

In the next two sections I present a brief background to ozone chemistry 
and describe current issues in stratospheric ozone. The last section of this Chapter 
outlines the contribution of this work to research ongoing in the scientific community. 

1.1. Background on Ozone Chemistry 

Ozone in the stratosphere is produced via photolysis of molecular oxygen: 

O2 + hi/ (A < 242 nm) —* 0+0 (1.1) 

0 + O2 + M — * O3 + M (1.2) 

where M is any molecule. Production is greatest over the tropics and low latitudes 
where the flux of ultraviolet radiation is greatest. Reaction 1.2 proceeds rapidly, as 
does the photolysis of ozone, which regenerates the oxygen atom: 

i 

O3 + hu — ► 0 + O2. (1-3) 

Throughout the daylit stratosphere, ozone is in fast photochemical equilibrium with 
atomic oxygen. The sum of ozone and atomic oxgen is called odd oxygen. A fourth 
reaction acts as a sink of odd oxygen: 

— * 2 O 2 . 


O3 + 0 


(1.4) 
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Reactions 1.1 through 1.4 are known as the Chapman reactions and represent a first 
attempt to describe the ozone budget (Chapman, 1930). 

The average motion of the stratosphere can be described as a kind of 
fountain, with warm air upwelling through the tropical tropopause, spreading out 
over the mid-latitudes, and then cooling and subsiding at high latitudes. This 
motion is known as the Brewer-Dobson circulation. The air above the tropical 
tropopause is photochemically “young” air — that is, air that has not been exposed 
long to the ultraviolet radiation levels typical of the stratosphere. Such air contains 
high levels of tropospheric gases — e.g., CH 4 , N2O, and the chlorofluorocarbons 
(CFCs) — and low levels of some of the radicals derived from these source gases 
— e.g., NO or CIO. As the air is pushed to higher altitudes and out over the mid- 
latitudes, the period of exposure to ultraviolet light lengthens, the source gases break 
down, and the abundances of radicals builds up. These radicals play important roles 
in catalytic ozone loss processes. 

One such cycle, involving odd hydrogen (OH-f HO 2 ), proceeds as follows: 

OH -f- O3 — ♦ HO2 ■+■ O2 (1.5) 

HO2 -f" O3 — ► OH -f 2O2 (1.6) 

net : 2 O3 — * 3O2. 

An important cycle involving chlorine and bromine radicals is: 


H02 + Z0 

— HOZ + 02 

(1.7) 

HOZ + hv 

- OH + Z 

( 1 . 8 ) 

Z + O 3 

— ► ZO O 2 

(1.9) 
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OH + O 3 — ► H0 2 + Oj (i.10) 

net : 2 O 3 — ► 3 0 2 


where Z represents either Cl or Br. Other halogen cycles include: 


BrO + CIO 
Br + O 3 
CI + O 3 
net : 2 O 3 

and 


Br + Cl + 0 2 

( 1 . 11 ) 

BrO + 0 2 

( 1 . 12 ) 

CIO + o 2 

(1.13) 

30 2 



CIO + O Cl + 0 2 (1.14) 

Cl 4 - O3 — + CIO ■}■ Oj (1.15) 

net : 0 -f* O 3 — ► 0 2 . 

Two cycles involving the nitrogen radicals NO and N0 2 are: 


NO + O 3 — ► NO 2 “h O 2 

(1.16) 

NO 2 *f 0 — ► NO + O 2 

(1.17) 

net : 0 3 + 0 — ► 2O2. 

* 


NO 4 * O 3 — ► NO 2 4* O 2 

(1.18) 

NO 2 + O 3 — ► NO 3 + O 2 

(1.19) 

NO 3 + hu — ► NO * 4 “ O 2 

( 1 . 20 ) 

net : 2 O 3 — > 3 O 2 . 
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Together the two nitrogen radicals are known sis NOx, or odd nitrogen. 

The task of identifying these cycles and determining their relative impor- 
tance has occupied atmospheric chemists since the 1960s, when calculations showed 
that the Chapman cycle led to an overestimate of ozone in the stratosphere. The 
nitrogen cycle, Reaction 1.16-1.17, was first described in the early 1970s (Crutzen, 
1970; Crutzen, 1971); the role of chlorine in the ozone budget was suggested shortly 
thereafter (Stolarski and Cicerone, 1974). More recently, measurements of radicals 
made from instruments onboard ER-2 aircraft have permitted careful analysis of 
the loss cycles in the lower stratosphere (e.g., Wennberg et al., 1994). The ER-2 
data show that in this region, the hydrogen cycle (Reactions 1.5- 1.6) and the halo- 
gen cycles (Reactions 1.7-1.10 and Reactions 1.11-1.13) dominate all other cycles. 
In the middle stratosphere, satellite data together with model calculations demon- 
strate the importance of the NOx cycles (Garcia and Solomon, 1994). Analysis of 
satellite data from the upper stratosphere show that the odd hydrogen cycle regains 
importance here (Jackman et al., 1986; Crutzen et al., 1995). 

Complicating the ozone budget are reactions that take place among the 
radicals themselves. For example, NO 2 can slow the chlorine catalytic cycles by 
tying up the chlorine radical CIO in the longer-lived, less reactive molecule CIONO2: 

CIO + N0 2 + M -> CIONO 2 + M. (1.21) 

The radical NO, on the other hand, influences the odd hydrogen cycle via this 
reaction: 


NO + H0 2 


NO 2 + OH. 


( 1 . 22 ) 
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This reaction competes with Reaction 1.6, the rate-limiting step in the odd hydrogen 
cycle. In other words, the two odd nitrogen species, NO and N0 2 , can short-circuit 
the catalytic cycles involving chlorine and odd hydrogen, leading to slower ozone 
loss rates. 

In the lower and middle stratosphere under most conditions, the produc- 
tion and loss rates of ozone are slow, on timescales of months. The long chemical 
lifetime of ozone means that its distribution is governed in large part by transport. 
Distinguishing between the effects of transport and chemistry on the ozone budget 
in these regions of the atmosphere is tricky, as is discussed in Chapter 5. In the 
upper stratosphere, above about 40 km, the flux of ultraviolet sunlight increases and 
production of ozone via the photolysis of molecular oxygen proceeds rapidly. Loss 
processes are also rapid due to the presence here of large abundances of the necessary 
radicals. In this region, ozone is considered to be in photochemical steady-state. 

1.2. Issues in Stratospheric Ozone 

Not long after it was first suggested that chlorine played an important role 
in the ozone budget (Stolarski and Cicerone, 1974), Molina and Rowland (1974) 
perceived the possible consequences that increasing emissions of CFCs could have 
on the stratospheric ozone layer. Molina and Rowland pointed out that the CFCs, 
though stable in the troposphere where the flux of ultraviolet sunlight is greatly 
attenuated, would photodissociate much more readily in the harsher radiation en- 
vironment of the stratosphere. Breakdown of the CFCs would increase abundances 
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of atomic chlorine and CIO radicals, which could then participate in cycles such as 
Reactions 1.14-1.15, accelerating ozone loss rates. 

In response to the ideas of Molina and Rowland, missions such as the 
Upper Atmosphere Research Satellite (UARS) were planned. These satellites were 
designed to monitor chemistry in the mid- to upper stratosphere over the tropics 
and mid-latitudes. Ironically, the first signs of ozone perturbation were reported 
by Farman (1985) over Antarctica. His discovery marked the beginning of intense 
research activity in polar regions as well as in laboratories worldwide. (See, for 
example, Solomon et al., 1986; McElroy et al., 1986; Toon et al., 1986; Crutzen 
and Arnold, 1986; Tolbert et al., 1987; and the special issues of the Journal of 
Geophysical Research , 94 , nos. D9 and D14.) By the early 1990s, an explanation 
of the annual loss of springtime ozone over Antarctica had been pieced together 
(Solomon, 1990). In brief, the theory states during the long, dark polar winter, 
reactions taking place on the surfaces of polar stratospheric clouds shift atmospheric 
chlorine from relatively unreactive forms to much less stable forms. At the same 
time, heterogeneous reactions on the cloud particles shift atmospheric nitrogen in 
the opposite direction, from reactive to reservoir species. The reservoir nitrogen, 
HNO3, may remain dissolved in the cloud particle. Over the Antarctic, where winter 

If 

temperatures drop to extremely low values and the cloud particles grow to large 
sizes, some reactive nitrogen may be carried permanently out of the stratosphere by 
particle sedimentation. In spring, when sunlight returns to the region, the less stable 
chlorine molecules quickly photolyze, leading to abundances of chlorine radicals 
many times that ususually encountered in the stratosphere. The low levels of reactive 
nitrogen mean that Reaction 1.21 proceeds slowly, and the abundances of reactive 
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chlorine may persist for periods of days to weeks. The reactive chlorine species 
include atomic chlorine and CIO, which together drive rapid catalytic cycles such 
as those described in the previous section, destroying ozone at the rate of about one 
per cent per day (Manney et al., 1993). 

Work continues on many issues concerning the annual loss of ozone over 
Antarctica. For example, the tie between stratospheric chlorine abundances and 
Antarctic ozone loss rates was strengthened by recent data from aircraft and satel- 
lites, thus implicating chlorine in ozone hole chemistry (Waters et al., 1993). Ozone 
loss during Arctic spring has also been carefully monitored (Webster et al., 1993; 
Profitt et al., 1993; Manney et al., 1994a). The composition of the particles that 
make up polar stratospheric clouds is under investigation, as are the uptake coef- 
ficients for gases that react in or on the surfaces of these particles (e.g., Abbatt 
and Molina, 1992; Molina et al., 1993; Middlebrook et al., 1993). One important 
issue to be settled is how readily air flows through the Arctic polar vortex in late 
winter and early spring, when the winds that set up this vortex act as a barrier 
against large-scale mixing with mid-latitude air (Tuck et al., 1992; Plumb et al., 
1994). A porous barrier could lead to the propagation of ozone-depleted air to lower 
latitudes. Over the Antarctic, the current consensus seems to be that mixing with 
mid-latitude air occurs only at tlie vortex edge, and that air deep within the vortex 
remains isolated until long into spring (Manney et al., 1994b; WMO 1995). Lastly, 
the chemical mechanisms that shut down the the ozone loss processes in spring are 
under review. 

Work also continues on analysis of the chemistry of ozone over the tropics 
and mid-latitudes. As mentioned above, the main focus of this research has been on 
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identifying the various ozone loss cycles, and determining their relative importance 
at different times of the year and in different regions of the stratosphere. Much 
progress in understanding ozone chemistry was made in the late 1980s, when the 
importance of heterogeneous chemistry occurring on the surfaces of sulfate particles 
was understood (Hofmann and Solomon, 1989b). Like particles in polar strato- 
spheric clouds, these particles provide a means to convert reactive nitrogen to the 
reservoir species HNO 3 , thus affecting not only the NOx catalytic cycles, but also, 
indirectly, the chlorine and odd hydrogen cycles. Sulfate particles are ubiquitous in 
the lower stratosphere, but their abundances can grow by many orders of magnitude 
in the aftermath of a volcanic eruption forceful enough to inject clouds of SO 2 into 
the stratosphere. The indirect role these particles play in the ozone budget has been 
the subject of much research, especially in the years following the eruption of Mount 
Pinatubo in 1991 (Fahey et al., 1993; Kawa et al., 1993; Mills et al., 1993). 

Much current research also focusses on assigning causes to the long-term 
decadal trends in ozone abundances observed over the mid-latitudes. As a first step, 
knowledge of the ozone budget has been tested in the middle and upper stratosphere, 
where ozone is thought to be in photochemical equilibrium. Since photochemistry, 
not transport, determines the the ozone distribution in this region, testing the bud- 
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get here is straightforward. Early satellite observations from the late 1970s and 
mid-1980s showed larger abundances of ozone than predicted (Jackman et al., 1986; 
Natarajan and Callis, 1989). In the 1990s, new ozone observations made by UARS 
provided an opportunity to revisit the issue of the ozone budget. While the calcu- 
lations of Crutzen et al. (1995) showed a good match with observations, those of 
Dessler et al. (1996) showed some discrepancies. Reasons for the different outcomes 
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of the two studies are not clear. It is not surprising then that much remains uncer- 
tain about the causes of the long-term ozone trends. Calculations constrained by 
plausible scenarios of anthropogenic chlorine emissions do not simulate the observed 
decline in mid-latitude ozone (Chandra et al., 1995; C. Jackman, personal commu- 
nication). These calculations include simulations of heterogeneous chemistry taking 
place on sulfate aerosol. While chlorine is certainly implicated in the ozone loss, 
other possible causes are under investigation. For example, the role of bromine and 
the possible role of iodine in the stratospheric ozone budget are both being studied 
(Garcia and Solomon, 1994; Solomon et ah, 1994; Avallone et ah, 1995; Wennberg 
et ah, 1996). 

Still another major issue involving stratospheric ozone focusses on the pos- 
sible impact that supersonic transports (also called high-speed civilian transports, 
or HSCTs) could have on stratospheric ozone abundances. The emissions of these 
aircraft include large quantities of reactive nitrogen oxides (NOy), water vapor, and 
SO 2 (which may be rapidly converted to sulfate aerosol). The HSCT issue, first in- 
vestigated in the early 1970s (Johnston, 1971), is still the subject of much research. 
Over the years, atmospheric chemists have debated the magnitude and even the sign 
of the change in stratospheric ozone abundances that would follow the deployment 
of a fleet of such aircraft. One problem that chemists encounter in their attempts to 
resolve this issue is that transport processes in the stratosphere are not well under- 
stood, and so the ultimate fate of the HSCT emissions is not known. For example, 
if mid-latitude air, rich in HSCT NOy emissions, can mix with tropical air, then the 
upwelling motions typical of the tropics may carry the NOy-perturbed air up to the 
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middle stratosphere (Avallone and Prather, 1996). In this region, as noted above, 
the odd nitrogen catalytic cycles are thought to dominate ozone loss processes. 

A last ozone issue concerns the temperature changes and radiative forcing 
effects that may accompany loss of stratospheric ozone or increases in greenhouse 
gases such as CO 2 . While the net effect of stratospheric ozone loss on the earth’s 
surface temperature is under debate, there is consensus that a thinner ozone layer 
would cool the stratosphere as would an increase of greenhouse gases (Ramaswamy 
et al., 1992; Hauglustaine et al., 1994). The chemical consequences of a colder 
stratosphere are currently being examined. For example, Pitari and Ricciardulli 
(1994) have speculated that colder temperatures over the Arctic in winter could 
lead to larger ozone losses there. 

1.3. This Work 

This work focusses on chemical processes important to ozone chemisty 
in the mid- and lower stratosphere. The research analyzes observations made by 
the Halogen Occultation Experiment (HALOE), an instrument onboard the UARS 
satellite. HALOE measures the abundances of a suite of gases important in ozone 
chemistry, including ozone itself. To solve the problem of distinguishing between 

•t 

chemical change and change due to transport, we make use of two tracers measured 
by HALOE. In brief, we use tracer analysis to “tag” each air parcel observed with an 
identity related to its photochemical age. To our knowledge, no previous work has 
followed the evolution of short-lived species over time with this system of tagging. 
We apply this technique of tracer analysis to two problems: first, we examine the 
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response of odd nitrogen and ozone to the 1991 Mount Pinatubo eruption, which 
injected huge quantities of sulfate aerosol into the stratosphere. This work confirms 
earlier predictions of the sensitivity of ozone to changes in NOx in the middle strato- 
sphere. Second, we examine the recovery phase of the Antarctic ozone hole during 
spring. The second work shows how reactive chlorine may be quickly converted to 
the reservoir species HC1 in the aftermath of severe ozone depletion, putting an end 
to further ozone loss. 


9 



t 


CHAPTER 2 
THE DATA SET 

Planning for HALOE began in 1975, shortly after Rowland and Molina 
(1974) first suggested that increasing emissions of CFCs could harm the stratospheric 
ozone layer. In the early 1980’s, HALOE was chosen as one of ten instruments to be 
included on UARS. The satellite was finally launched in September 1991, and the 
instrument began making measurements that October (Russell et al., 1993a). As of 
this writing, HALOE continues observations. 

HALOE takes advantage of the geometry of occultation to look through 
the atmosphere toward the sun at local sunrise and sunset. (See Figure 1.) The 
instrument determines the abundances of seven gases — ozone, NO, NO 2 , HC1, HF, 
CH 4 , and H 2 O — by measuring the absorption of infrared sunlight within specified 
wavelength intervals. These seven gases were carefully chosen. Measurements of 
ozone provide information on long- and short-term trends of this important molecule. 
Observations of NO and NO 2 are useful since these radicals play direct and indirect 
roles in the ozone budget. Observations of CH 4 together with those of H 2 O can be 
used to derive abundances of OH, another radical important in the ozone budget. 
The molecule HC1 was chosen as part of the observation set since it is a reservoir of 
atmospheric chlorine. Since all HF in the stratosphere derives from CFCs, HALOE 
measurements of HF can provide information on the influx of CFCs up through 
the tropopause. Tracking HC1 and HF together can tell us how much the chlorine 
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content of the stratosphere has increased due to increasing emissions of the CFCs. 
In addition, HF and CH4, both very long-lived molecules with lifetimes on the order 
of years through much of the range of HALOE, serve as useful tracers of atmospheric 
motion. 

Concentrations of an eighth molecule, CO2, are also measured. As is de- 
scribed below, the measurements of C0 2 provide a means to determine the altitude 
grid of the profiles, as well as atmospheric temperature and pressure. 

In the following sections, I discuss in further detail the following topics: 
(1) observation geometry and geographic coverage, (2) measurement techniques, (3) 
determination of pressure and temperature, (4) sources of error, and (5) efforts to 
validate the measurements. 

2.1. Observation Geometry and Geographic 

Coverage 

UARS circles the globe 15 times every 24 hours in an orbit inclined 57 
degrees to the equatorial plane. During occultation, when the atmosphere comes 
between the satellite and the sun, HALOE’s telescope aims toward the sun and 
obtains a vertical scan of the atmosphere. The time of the measurement coincides 
with sunrise or sunset at the tangent point on the earth’s surface. The instrument 
makes about 30 scans a day, corresponding to 15 sunrise and 15 sunset events. Fig- 
ure 1 shows the geometry of the solar occultation experiment. HALOE observations 
range from the lower stratosphere up to about 150 km. For those species with 
insignificant concentrations in the mesosphere and thermosphere, HALOE profiles 
reach an upper limit of about 50 km. Below the tropopause, water vapor and clouds 
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cause considerable interference in the measurements. The vertical resolution of the 
instrument ranges from 2 to 4 km, depending on the channel, and the horizontal 
resolution is about 300 km along the limb. 

Since the earth is turning slowly during the course of a day’s measurements, 
the sunrise measurements for each day fall at different locations along a narrow 
latitude band, as do the sunsets. In addition, the orbit of the satellite precesses 
slowly along the earth-sun line, with one complete revolution around the earth-sun 
line taking about 8 weeks. This precession enables the instrument to sample sunrises 
and sunsets within a different latitude band each day. Precession of the satellite’s 
orbit together with the earth’s rotation and revolution about the sun result in a 
complicated measurement pattern that spirals up and down the earth. Figure 2 
shows the measurement pattern for the first 30 days of operation, from October 10, 
1991, to November 11 of that year. The circles represent sunrises and the crosses, 
sunsets. As can be seen from the plot, the sunrise and sunset measurements during 
this time period together covered most of the globe, from about 55 N to 80 S. It 
is also apparent in the plot that the measurement points tended to cluster at high 
latitudes. 

To derive vertical profiles of the species measured at local sunrise and 
sunset, an “onion peel” method is employed. The method takes advantage of the 
curvature of the atmosphere to build a profile at the tangent point. To understand 
this method, it is helpful to imagine the atmosphere divided up into very thin shells 
or layers. The onion-peel method rests on the assumption that the concentration 
of each gas measured is uniform across each thin shell. At the top of a HALOE 
scan, the viewing limb crosses through shells containing little or no target gas. As 
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the instruments scans deeper into the atmosphere, the viewing limb finally reaches 
the uppermost shell of the gas — that is, the first shell that contains measurable 
quantities of that gas. At that time, nearly all of the viewing limb crosses shells 
in which where the species is not present, and the absorption signal is due entirely 
to the target gas located at the tangent point. A calculation to determine the 
concentration of the target gas at the tangent point may then be performed. A 
little lower in the atmosphere, the limb passes through the empty shells, part of the 
uppermost shell, and a small increment of a new shell at the tangent point. Since 
the target gas concentration of the uppermost shell is now known, the distribution 
of the gas along the limb can be determined together with the concentration of the 
new shell at the tangent point. In this way, working increment by increment from 
higher altitudes to lower ones and constructing distributions of gases along each 
viewing limb, vertical profiles of each gas can be obtained. 

The solar occultation technique has several advantages over nadir-viewing 
instruments that may depend on, for example, radiative emission of molecules. First, 
the abundance of the absorbers along the limb path is 30 to 60 times that along the 
viewing path of a downward-looking instrument. The thick optical depth enables 
the instrument to measure with greater accuracy molecules in small abundance such 

ST 

as NO or NO 2 . Second, the occultation geometry also allows for higher vertical res- 
olution. Third, calibration of the instrument can easily be made by pointing the in- 
strument toward the sun during that part of the satellite’s orbit when no atmosphere 
intervenes between the instrument and the sun. Lastly, the solar occulation tech- 
nique requires somewhat less precise knowledge of atmospheric temperature than do 
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techniques relying on radiative emission. Radiative equilibrium in the atmosphere 
is not a necessary condition for observation. 

The technique also has several disadvantages. To begin with, the assump- 
tion that each thin shell of the amosphere has a uniform concentration can be flawed. 
In regions of high horizontal gradients, the concentrations of gases across the shells 
may vary considerably. For most molecules under most conditions, however, errors 
arising from this assumption are minimal since the decreasing density of the atmo- 
sphere with height means that the contribution to the signal is usually greatest at 
or close to the tangent point. The uniform-shell assumption leads to greatest error 
for those molecules that have a diurnal cycle — NO and NO 2 in this case. During 
sunset, NO converts rapidly to N0 2 ; at sunrise, the opposite takes place. While 
no part of the observing limb lies in darkness, the solar zenith angle varies along 
the limb, giving rise to distributions of NO and N0 2 that are not symmetric about 
the tangent point. To take this asymmetry into account, the measurements of NO 
and NO 2 are adjusted using known chemical theory (Boughner et al., 1980). Above 
about 30 mb, the adjustment is small, resulting in changes in NO or N0 2 mixing ra- 
tios of less than 5 per cent. Below 30 mb, the adjustment is much greater, resulting 
in changes up to 50 per cent. The larger adjustment at lower altitudes is necessary 

c 

since in this region the viewing limb crosses through so many layers of atmosphere, 
including layers containing large abundances of NO and N0 2 . 

Another disadvantage of the solar occultation technique is that on any 
given day it provides information over only a limited part of the globe. Figure 3 
shows the daily average latitude of sunrise and sunset measurements during late 
spring and early summer of 1992. Again the circles represent sunrise observations, 
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and the crosses, sunset. As is clear from the plot, a picture of the atmosphere as a 
whole must be pieced together from a succession of measurements. This piecemeal 
approach to observation presents less of a problem when tracking short-term trends 
at high latitudes during those intervals when the instrument repeatedly samples the 
same region several days in a row. 

2.2. Measurement Techniques 

HALOE measures the abundances of trace gases in the atmosphere via 
two methods — the gas filter correlation method and broadband radiometry. Both 
methods involve measuring the absorption of infrared sunlight taking place in the 
atmosphere within certain narrow wavelength intervals. The wavelength intervals 
chosen are spectral regions where the target molecules absorb strongly. The central 
wavelengths of these intervals range from 2.452 microns for HF to 9.852 microns for 
ozone; bandwidths range from 15 cm" 1 for NO 2 to 100 cm -1 for ozone. With inde- 
pendent knowledge of the absorbtion spectra of these molecules, one can calculate 
the distribution of each species present along the limb. The gas filter correlation 
method is the more sensitive of the two methods, and is used to measure NO, HCl, 
HF, and CH 4 . Broadband radiometry is used for ozone, N0 2 , H 2 0, and C0 2 . Both 

t 

methods are discussed in greater detail below. 

HALOE measures the abundances of all eight traces gases simultaneously 
during sunrise and sunset events. Sunlight entering the instrument’s telescope is 
split by a chopper that is gold-plated and polished (and thus highly reflective) on 
one side. Half the solar energy is reflected off the chopper toward the broadband 
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radiometers, and half is directed toward the gas filter correlation channels. Figure 4 
shows the layout of the optical train. Beam splitters, marked BS in the diagram, 
route the sunlight to the different components. Broadband filters, labeled BF, cut 
out undesired wavelengths of light so that the energy reaching each of the detectors 
belongs to a narrow wavelength interval, between 44 and 110 cm -1 in width. The 
central wavelengths of the intervals measured by the eight channels are marked on 
the diagram. For example, the central wavelength for the NO channel is 5.263 
microns. 

2.2.1. Gas filter correlation technique 

The absorption spectra of many molecules in the atmosphere overlap to 
some degree, and isolating the absorption bands of one species can be difficult. 
For example, CH* is a major interfering species in the absorption spectra of HC1. 
The problem is especially acute since CH4 is 1000 times more abundant in the 
atmosphere than is HC1. Filtering out such interferences requires special techniques. 
It should be kept in mind that the detectors on HALOE “see” an entire wavelength 
interval at once; the instrument is not capable of scanning slowly through an interval, 
wavelength increment by increment. 

r 

To overcome the problem of interfering species, HALOE employs the gas 
correlation method. (See Figure 5a and b.) In this method, incoming sunlight in 
each channel is split into two beams. One beam passes through a reference cell 
containing the target gas and then on to a detector; the other beam passes through 
a vacuum and thence to the detector. Calibration takes place once each orbit, when 
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the satellite travels between the earth and the sun, and the instrument can look 
directly toward the sun, without any atmosphere intervening. At this time, the 
detector at the end of the gas-cell path records the transmission of sunlight through 
the cell: 


Tctii = e- T » (2.1) 

where T g is the effective optical depth of the cell integrated over the wavelength 
interval. At the same time, transmission through the vacuum path, T vacuum , is, of 
course, unity. A gain factor G is then applied to the signal arising from the gas-cell 
path so that transmission through the two paths balance: 


GTcell =S Tvacuum = 1. (2.2) 

During occultation, when the instrument encounters the atmosphere in its 
view toward the sun, the transmission of sunlight along the two paths changes. The 
detector at the end of the vacuum path sees the attenuation of infrared sunlight due 
to the target gas within the atmosphere as well as that due to interfering gases: 


T 

* vacuum 


= e -T ° e~ T \ 


(2.3) 


where r a is the optical depth of the species of interest and r,- is the optical depth 
of the interfering species. The detector at the end of the gas-cell path, on the 
other hand, sees a different transmission. Since the sunlight that travels this path is 
already strongly attenuated by the gas enclosed inside the reference cell, the detector 
cannot, in fact, see any additional absorption taking place due to the target species 
within the atmosphere. (Absorption of sunlight is linear with absorber abundance 
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only in the limit of very small concentrations of absorbing material.) Within the cell, 
absorption processes become saturated — that is, the incoming sunlight excites the 
gas cell molecules to higher vibrational levels, but the relaxation processes taking 
place cannot keep up with the continued influx of energy. The detector in the 
gas-cell path can, however, monitor the attenuation of sunlight due to interfering 
species since absorption within these bands is not saturated. Transmission through 
the gas-cell path may therefore be written: 


T ce ii = G e Tg 


( 2 . 4 ) 


Abundance of the species of interest along the limb may be calculated 
by taking the ratio of the difference in transmission along the two paths to the 
transmission along the vacuum path: 


T ce u — T v , 


T 

x vacuum 


G e -T » e~ Ti - e -T “ e“ T ' 


9 — j9'“l’1 


( 2 . 5 ) 


= 1 - e~ Ta 


( 2 . 6 ) 


= j _ e -/p(0 ff (0 dl 


( 2 . 7 ) 


Here p is the density of the absorbing species, a is the cross section of the species, 
and dl is a horizontal increment along the limb. The gas cell in this way acts as 

C 

filter to remove the effects of the interfering gases within the wavelength interval, 
even when absorption due to the interfering gas is many times greater than that due 
to the species of interest. 

A slightly different set-up is used to measure HC1. Since CH 4 absorbs so 
strongly in the HC1 wavelength region, an extra cell containing CH 4 intercepts the 
sunlight before it is split into the HC1 gas-cell path and the vacuum path (Figure 5 b). 
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In this way, a significant part of e~ T ' is pulled out of the subtraction in the numerator 
on the righthand side of Equation 2.6, and the ratio on the lefthand side of the 
equation is calculated with greater accuracy. 

2.2.2. Broadband radiometry 

Measurement of relatively abundant species in the atmosphere can be made 
more simply than the above method. Both ozone and H 2 O are present in large 
enough quantities and absorb so strongly within their specified wavelength intervals 
that under most conditions we can assume that the absorption taking place within 
the intervals is due entirely to ozone or to H 2 O. In other words, the absorption due to 
interfering gases in these channels is minimal, and a simpler method of measurement, 
broadband radiometry, may be employed. (See Figure 5c.) A third gas measured 
by broadband radiometry, NO 2 , is not in fact present in large quantities in the 
atmosphere, but under normal conditions this gas is the dominant absorber in the 
spectral region chosen. Only water as a minor but significant interferent. 

In broadband radiometry, the detectors are first calibrated by measuring 
the intensity of sunlight during that time of the orbit when no atmosphere lies 
between the satellite and the sun. The ratio of that intensity, Jo ? to the intensity of 

«r 

light attenuated by the atmosphere during occultation gives the abundance of the 
absorbing gas along the limb: 





( 2 . 8 ) 

(2.9) 



24 


Since H 2 0 and CH 4 absorb in the same wavelength region as N0 2 , the abundance 
of these interfering gases along the limb are first calculated and the effect of these 
abundances on the N0 2 signal determined. In this way, the absorption due solely 
to N0 2 may be evaluated. 

2.2.3. Measuring pressure and temperature 

HALOE measurements of C0 2 provide a means to determine the pressure 
and temperature profiles for each sunrise and sunset event, and to assign an altitude 
to each data point. The algorithms developed to calculate these variables take 
advantage of the fact that C0 2 is extremely well-mixed througout much of the 
atmosphere, with a nearly constant mixing ratio of 350 ppm. 

As a first step, the altitude of each data point in a profile is estimated 
from flight path information together with the apparent oblateness of the sun at 
the time of measurement. To determine the altitude with greater certainty, we 
next make use of measurements supplied by the National Meteorological Center 
(NMC). In the 29-39 km region, the measured absorbance of sunlight by C0 2 at 2.8 
microns at each estimated altitude gridpoint is compared to the C0 2 absorbance 

calculated from NMC pressure and temperature profiles. The 29-39 km region was 

* 

chosen as one in which both HALOE and NMC make measurements with reasonable 
accuracy. Differences between the observed and calculated absorption are minimized 
by adjusting the estimated altitude grid of the HALOE observations. Once the two 
CO 2 profiles match, the altitude grid of the HALOE data between 29 and 39 km is 
known, as well as the “sink rate” of HALOE observations — that is, the distance up 
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or down in the atmosphere the HALOE observations shifts per increment of time. 
Since the sink rate is nearly constant during an event, the altitude of every data 
point may be easily determined by extrapolation from the 29-39 km region. 

Next, knowledge of the absorption due to COj at 2.8 microns allows us to 
calculate the concentration of CO 2 along the tangent profile. Since the mixing ratio 
of CO 2 is nearly constant up to about 85 km, a profile of total air density is readily 
calculated up to that altitude. Some adjustment of the total air density is made 
to take into account the increasing abundance of CO 2 due to increasing emissions 
of hydrocarbons at the earth’s surface. Lastly, profiles of pressure and temperature 
up to 85 km can be constructed using the hydrostatic approximation and the ideal 
gas law. Above 85 km, a model based on climatology is used to determine these 
quantities. 

2.3. Sources of Error 

Sources of random error in the HALOE observations of trace species in- 
clude instrument noise, sun-tracker jitter, and atmospheric temperature uncertain- 
ties. Uncertainty in temperature leads to uncertainty in interpreting the absorption 
spectra. Random error may also propagate from one channel to another when calcu- 

c 

lating the effects of interfering gases. Sources of systematic error include uncertainty 
in absorption line strengths and the reference pressures calculated with CO 2 . Other 
sources of systematic error include the error introduced when calculating the vari- 
ation of the diurnally changing species NO and N0 2 along the viewing limb. In 
regions of steep horizontal gradients, systematic error may also be caused by as- 
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suming that the distributions of the longer lived species are symmetric along the 
viewing limb. 

To evaluate these errors, workers at NASA Langley performed a Monte 
Carlo analysis which repeatedly simulated the retrieval of mixing ratios, introducing 
a new set of errors with each iteration. These calculations predict 10 to 15 per 
cent error in all species between about 30 and 1 mb. Above 1 mb, concentrations 
fall rapidly with height, increasing signal-to-noise ratios and thus error. Below 30 
mb, retrieval of the mixing ratios is made difficult by the increasing opacity of the 
atmosphere along the limb and, for some species, small concentrations typical of 
the lower stratosphere. Under normal conditions, per cent error for most molecules 
above 1 mb and below 30 mb ranges from 20 to 100 per cent. During the early years 
of the mission, however, observations at the lower altitudes were made especially 
uncertain by the presence in this region of large abundances of sulfate aerosol leftover 
from the June 1991 eruption of Mount Pinatubo. Under heavy aerosol conditions, 
relative error could increase to several hundred per cent for NO 2 , measured by the 
broadband radiometry method. Meausrements of the other broadband-radiometer 
species were also affected, though not to the same degree. Using information from 

the gas-correlation channels, transmission in the NO 2 channel was corrected to take 
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aerosol absorption into account. 


2.4. Validation Efforts 

Efforts were made to validate the HALOE dataset — that is, we attempted 
to judge the quality of the data by comparing HALOE data to measurements made 
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by other instruments and by performing consistency checks. Our group focussed on 
validating observations of HALOE NO and NOj. 

2-4- 1- Comparison to other measurements 

I present here a brief overview of the work done to compare HALOE NO 
and NO 2 to observations of these molecules made by other intstruments. The reader 
is referred to Gordley et al. (1996) for a more detailed description of this effort. 

Because of the short chemical lifetimes of NO and NO 2 (on the order of 
seconds), stringent requirements were made on the time and space coincidence of 
pairs of observations. For example, Figure 6 compares a HALOE N0 2 profile to 
a profile measured by a balloon-borne interferometer (Toon, 1991). The HALOE 
measurements were made September 25, 1993, and the balloon measurements, the 
next day. The two profiles are located within 1 degree latitude and 0.5 degrees 
longitude of each other, at about 34 N. Observations made by the interferometer 
are represented by diamonds; HALOE measurements are denoted by stars. The 
two profiles show relatively good agreement, though HALOE observations fall short 
of the balloon observations at the maximum and below about 40 mb. A similar 

comparison for NO is shown in Figure 7. Again the two profiles match fairly well. 

€ 

As a second example, Figure 8 compares sunrise HALOE NO 2 total col- 
umn amount to measurements made by a ground-based spectrometer at Fritz Peak 
Observatory in Colorado for the years 1992 and 1993. Figure 9 shows a similar com- 
parison for sunset data. In both plots, the circles and triangles represent HALOE 
data for 1992 and 1993, respectively. The stars and dashes represent the two years 
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of ground-based measurements. All of the HALOE data included in the two plots 
fall within 6 degrees latitude and 12 degrees longitude of Fritz Peak, located at 40 N 
and 106 W. The plots show the seasonal maximum for both years occurs in summer, 
when HNO3 photolysis is fastest, and more N0 2 forms. For both sunset and sunrise 
measurements, HALOE tends to show low amounts of column NO2 relative to the 
Fritz Peak measurements. 

Analyzing several such comparisons for NO and NO2, we came to the 
following conclusions. First, both HALOE observations of both NO and NO2 show 
relatively good agreement with observations made by other instruments. HALOE 
NO2 in the lower stratosphere, below about 40 mb, tends to be somewhat low, 
perhaps due in large part to an overcorrection for sulfate aerosol absorption. (See 
Section 2.3). Another possible source of error at low altitudes may be uncertainties 
in the CO2 retrieval in a region where significant CO2 line broadening occurs. Errors 
in interpreting the CO2 signal would lead to errors in calculating total air density 
and pressure. For NO2, whose mixing ratio changes rapidly with height in the 
lower stratosphere, accurate knowledge of air density and pressure is critical in the 
retrieval. It is conceivable that an underestimate of NO2 in the lowest few kilometers 

of the stratosphere may have led to the low NO2 column amounts noted in Figures 8 

» 

and 9. 


2.^.2. Sunrise to sunset comparisons 


An important internal consistency check on HALOE NO and NO2 data 
may be done by comparing sunrise and sunset profiles at times in the mission when 
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the two narrow latitude bands of sunrise and sunset observations overlap. This 
coincidence of sunrise and sunset measurements occurs about 6 to 7 times a year, 
usually over low latitudes. Comparing measurements of NO and NO2 at sunrise 
and sunset requires some care since the concentrations of both species are changing 
rapidly at these times. During the day, a set of fast chemical reactions establishes 
photochemical equilibrium between the two species. But at sunset, the reactions 
that during the day convert NO2 back to NO — i.e., the photolysis of NO2 and the 
reaction of NO2 with atomic oxygen — slow down and finally cease. NO concentra- 
tions quickly fall by reaction with ozone, and NO2 quickly rises. At sunrise, as NO2 
once again begins to photolyze and the concentrations of atomic oxygen increase, 

NO abundances rapidly recover and NO2 declines. Because of the fast chemistry 

r 

occuring between NO and NO2 at sunrise and sunset, it is best for this task to look 
at the sum of the two species, NO+NO2. 

Still another factor to consider in making sunrise to sunset comparisons is 
that during the night NO2 is slowly converted to N2O5 through the reactions: 

NO2 "f O3 — * NO3 -|- O2 (2.10) 

NO2 +N 0 3 +M — N 2 O s +M (2.11) 

m 

where M is any molecule. The rate-limiting step in this pair of reactions is Re- 
action 2.10, and NO3 is considered to be in steady state. The rate of decay NO2 
during the night may then be written: 


d ln[N0 2 ]/dt = — 2Jfc[0 3 ] 


( 2 . 12 ) 
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where the reaction rate k is a function of temperature, T (DeMore et al. 1992), 

k = 1.2 x 10 _13 exp(— 2450/7") cmV 1 . (2.13) 

Using HALOE O3 data and temperatures, one may calculate the expected sunrise 
NO+NOj values from the measured N0+N0 2 at sunset. The predicted sunrise 
concentrations can then be compared to the sunrise sums of NO+NOj as measured 
by HALOE. 

In carrying out these comparisons, pairs of sunrise and sunset events were 
chosen such that the sunrise event took place the morning following the sunset event 
in a location not more than 3 degrees latitude and 10 degrees longititude distant from 
the sunset. Ozone values were taken from the sunset events as were the temperature 
profiles. Figure 10 shows a sample comparison of HALOE sunrise to model sunrise. 
The data are from November 9, 1993, at about 17 N. It can be seen in Figure 10 
that the match between predicted and measured NO+N0 2 at sunrise is quite good. 
This agreement shows the internal consistency of HALOE’s NO and N0 2 , as well 
as the instrument’s ozone measurements. 

The above calculation was performed for many pairs of sunrise and sunset 
profiles with similarly good results, especially over low latitudes. For the few pairs 
at mid-latitudes, agreement was# not as good for reasons that are unclear. Above 
about 1 mb, HALOE measures little or no N0 2 at sunrise or sunset. In the meso- 
sphere at night, NO converts to N0 2 but the N0 2 photolyzes extremely rapidly once 
daylight arrives. Insufficient ozone is available at these altitudes to convert signif- 
icant amounts of nighttime N0 2 to N 2 0$. Therefore it comes as a surprise that 
HALOE tends to measure low NO at sunrise as compared to sunset. This mismatch 
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of sunrise to sunset at the higher altitudes may be due to some unknown, nighttime 
chemical process or perhaps simply to problems measuring NO abundances at these 
altitudes where concentrations are very small. 


r 
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Figure 1 . Solar occultation experiment geometry. 
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Figure 2. Latitude of HALOE observation for the first month of operation. The 
circles represent sunrise; the stars, sunset. 
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Figure 3. Average daily latitude for late spring and summer of 1992. The circles 
represent sunrise; the stars, sunset. 
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4 Gas Correlation Channels 
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Figure 4. HALOE optical diagram (Russell et al., 1993a). 
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ATMOSPHERE 
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Figure 5. HALOE measurement techniques, (a) Gas filter correlation technique, (b) 
Special gas filter correlation technique for HC1, (c) Broadband Radiometry (Russell 
et al., 1993a). 
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Comparison of Balloon Observations to HALOE 
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Figure 6. Comparison of balloon NO 2 measurements to HALOE NOj. The stars 
represent HALOE data, and the diamonds, balloon data. HALOE data are from 
Sept. 25, 1993, at 33 N and 110 W. Balloon data are from Sept. 26, 1993, at 34 N 
and 110 W. 
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Comparison of Balloon Observations to HALOE 



Figure 7. Same as preceding figure, but for NO. 
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Sunrise Column N02 in 1992-1993 



Figure 8. Ground-based observations of total column NOj at sunrise compared to 
HALOE total column NO 2 . Column abundances are in mol/cm 2 x 10 15 . The circles 
and triangles represent HALOE data from 1992 and 1993, respectively. The stars 
and equal signs denote the ground-based data for the two years. For location, see 
text. 
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Sunrise Column N02 in 1992-1993 
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Figure 9. Same as preceding figure, but for sunset. 
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NO + N02 Sunrise and Sunset Consistency Check 



Figure 10. Comparison of model to HALOE sunrise NO+NOj profiles. The model 
was initialized with HALOE sunset observations of ozone, NOj, and NO from the 
previous day. Data are from November 9, 1993, at about 17 N. 



CHAPTER 3 


MODELING STRATOSPHERIC CHEMISTRY 


The basic goal of most chemical modeling efforts is to predict the evolution 
of chemical species over time. The key equation to be solved at each timestep is the 
mass continuity equation: 


s -± = p-i - 8 J±- 8 J±- 6 Jjl 

St C Sx Sy Sz ’ 


( 3 - 1 ) 


where c is the concentration of a species, t is time, P is the production rate of c, L 
is the loss rate, and F x , F y , and F z are the fluxes of c in the zonal, meridional, and 
vertical directions, respectively. One major difficulty in solving such an equation 
is that it represents just one in a coupled system of partial differential equations, 
with each equation describing the time rate of change of a different species. The 
entire system must be solved simultaneously at each time step. Modeling transport 
processes together with chemistry also presents difficulties. In addition, in order to 
evaluate the photolysis rates of reactive species, the flux of solar radiation through 
the atmosphere must be properly calculated. 


For certain problems in the stratosphere, the mass continuity equation can 
be simplified somewhat and some or all of the flux terms dropped. For example, if 
our goal is to follow the chemistry occurring within an air parcel on the time scale of 
days, we may in some cases neglect mixing processes and utilize a simple box model 
describing the chemistry evolving at a discrete point. The point may be moved over 
time as the parcel is advected through space. Such an approach is taken for the 
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study described in Chapter 6 in which our concern is to track chemical change in a 
parcel of air within the Antarctic polar vortex in spring. The box-model approach is 
especially attractive when the goal of the study is to detect atmospheric change that 
is due solely to chemical processes, and not due to transport. For many studies that 
do require taking transport processes into account, a two-dimensional model of the 
atmosphere suffices. This is especially true in the stratosphere, where zonal mixing 
occurs rapidly, on timescales of days or weeks. The rapid zonal mixing ensures that 
the air along a latitude band is well-mixed and that the derivative of the flux with 
respect to x in Equation 3.1 may be dropped. We note that in regions such as the 
polar vortex where the horizontal gradients of species do not coincide with latitude 
circles, zonal averaging may lead to large error. 

For this work, we have employed two models — a box model and a two- 
dimensional model. The models are discussed in the next two sections. We describe 
the box model in fuller detail since it is used more extensively for this research. 

3.1. Box Model 

The box model — previously described by Logan et al. (1978), Prather 

and Jaffe (1990), and Salawitch et al. (1994) — contains a set of 35 chemical species 

«■ 

t 

and about 220 chemical reactions. The model employs a Newton-Raphson iterative 
method to solve implicitly for the mixing ratios of the species with short and inter- 
mediate lifetimes. Most of the reaction rate constants and absorption cross-sections 
in the model are taken from DeMore et al. (1992). Input parameters include a tem- 
perature profile, aerosol surface area density, latitude, and solar declination angle. 



The model has a variable timestep, ranging from tens of second at sunrise or sunset 
when many species are changing concentration rapidly to a couple hours during the 
night or at midday when chemical change is slow. The model may be initialized at 
either dawn or noon. 


3.1.1. Newton- Rafson method 


The Newton-Rafson method is useful for solving nonlinear systems of par- 
tial differential equations (Press et al., 1986). To understand this method, it is 
helpful to consider first just one equation in the system, i.e., 


0 = 


St 

= f(c ), 


+ P - Lc, 


(3.2) 

(3.3) 


which is simply Equation 3.1 rearranged and without the flux terms. The approach 
for solving this equation for c involves finding the slope of /(c) at some initial guess 
c,-, extrapolating the line given by the slope to the abscissa, and taking the value 
of c at the zero-crossing as the next guess, c,+i. The process is then repeated with 
the new c, + i. The Newton-Rafson method is analogous to finding the Taylor series 
about some point c: 

f(c + 6) ~ {(c) + 6 {'(c) + 6’Qfl + ... (3.4) 

where <5 is the correction to the current guess, and c + 6 is the next guess. For 
small values of <5, Equation 3.4 ~ 0, and we may neglect the higher order terms in 
Equation 3.4 so that 
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In this way, the correction S for each guess may be readily calculated. 

The advantage of the Newton-Rafson method is that the trial solutions 
converge rapidly onto the roots of the equations. It can be shown that having 
determined a correction 6{, the next correction, <5,+i, is 

(X6) 

For delta normalized to an initial guess, the squared dependence in this recurrence 
relation means that the guesses converge quadratically onto the solution. 

Taking the derivative of /(c) is straightforward. For example, consider a 
simplified continuity equation for NO: 

/([NO]) = 3^1 - /[NO,] + *[ 0 ,][NO] (3.7) 

where the brackets denote concentration, J is the photolysis rate of NO 2 , and k is 
the reaction rate constant of O 3 reacting with NO. The first term on the righthand 
side of the equation represents production of NO via NO 2 photolysis, and the second 
term, loss via NO + O 3 — * NO 2 + O 2 . For this continuity equation, the derivative 
of / with respect to [NO] can be written 

/'((NQJ) - + *[ 0 3 ][NO] (3.8) 

What complicates the calculation of photochemistry is that some 30 species 
in the model are changing with time, and the Newton-Rafson technique must be used 
for all these species simultaneously. In other words, roots of the partial differential 
equations must be found such that the solutions work for all species in the model. 
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Each function /, for each species c, may in fact be a function of many variables, 
such that 


f 1> 1 •■•Cn) 0* 


(3.9) 


The Taylor series expansion of /, in the neigborhood of the set of guessed concen- 
trations is: 


MC + SC) = MO + £ ^-Sc, + ... 

j=l 0C 3 


(3.10) 


where C is the vector of concentrations and SC is the vector of corrections. If the 
corrections are small, /*(C + SC) cz 0, and we may rewrite Equation 3.10 as 


where 


Aij6cj = -ft 
j=i 


Al] ~ Scj 


(3.11) 


(3.12) 


Equation 3.11 may be readily solved for the vector of corrections, 6C. These correc- 
tions are then added to the guessed vector C , and the process repeated with the new 
set of guesses. The iteration is halted when some degree of convergence is obtained. 

The model may be run 'in one of two modes. In the first mode, the model 
cycles through each model day for a specified time period and calculates the day-to- 
day evolution of all species over the time period. This approach is useful for tracking 
short-term chemical change, on the order of days or weeks. In the second mode, the 
model day is repeated over and until diet steady state is achieved — that is, until 
the the production and loss rates for each species integrated over the day balance 
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and the output gives reproducible diurnal cycles for all species. This approach 
is useful for determining the steady-state chemical response of the atmosphere to 
some long-term perturbation such as the increase in the sulfate-aerosol burden in 
the stratosphere that may follow an especially forceful volcanic eruption. 

3.1.2. Calculation of reaction rates 

As noted above, nearly all of the reaction rate constants and absorption 
cross-sections in the model are those recommended by DeMore et al. (1992). The 
absorption cross-sections for HNO 3 are those of Burkholder et al. (1993), and a 
quantum yield of 0.5 is used for both the CIO+NO 2 and the CI+NO 3 channels of 
CIONO 2 photolysis (Minton et al., 1992). The input temperature profile permits 
calculation of total air density, assuming hydrostatic equilibrium and estimating 
scale heights. The total air density profile may also be constructed from an input 
pressure profile together with temperature. Given the local temperature and density, 
all the reaction rates may be calculated. 

To calculate photolysis rates, first the flux of ultraviolet light through the 
atmosphere must be determined. The model is divided up into plane-parallel layers, 
and the flux through each layer^s computed by considering the direct and scattered 

t 

beam separately. A correction is made to take into account the sphericity of the 
earth’s atmosphere, and for twilight the flux of light scattered over the horizon is 
estimated for solar zenith angles up to 96 degrees (Prather, 1974). Two absorbers 
of ultraviolet radiation are considered: molecular oxygen and ozone. The profile of 
molecular oxygen is calculated from the total air density, and the profile of ozone 
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density is included as input. The temperature-dependent cross-sections of oxygen 
are those of Minschwaner et al. (1992 and 1993). 

Heterogenous chemistry in the model is carried out with rate constants 
given by the following (Cadle et al., 1975): 



where 7 represents the reaction efficiency (also known as the sticking coefficient), 
is the surface area of the particles per cubic centimeter of air, and v is the mean 
speed of the gas-phase reactant. To calculate the reaction rates of chlorine species on 
cold sulfate particles, the method developed by Hanson et al. (1994) is implemented. 

3.2. Two-Dimensional Model 

The two-dimensional model was developed by Paul Cruzen and co-workers 
at the Max Planck Institute for Air Chemistry in Mainz, Germany (Crutzen and 
Gidel, 1983; Briihl and Crutzen, 1988; Prather and Remsberg, 1993). The model 
calculates the chemical composition of the stratosphere and troposphere along con- 
stant pressure surfaces spaced about 2 km apart. It ranges in altitude from the 
earth’s surface to about 62 km and has a 10-degree latitude grid. The timestep 
of the model is usually set to be two hours long. The model includes about 120 
gas-phase and heterogeneous chemical reactions and about 50 photodissociations. 
Changes in chemical concentrations due to transport or chemistry are determined at 
each timestep. In order to determine the photolysis rates, a set of 12 radiation fields 
in the ultraviolet and visible range is recalculated once every 15 model days, one 
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field for each two-hour timestep in a day. For information on the model calculation 
of ultraviolet flux, see Zdunkowski (1980). 

Transport of chemical species in the model is split into two parts — trans- 
port by the net mean meridional and vertical winds and transport due to what is 
commonly called eddy diffusion. The net mean winds, which arise from temperature 
gradients and diabatic heating and cooling, are calculated off-line using temperature 
fields and heating rates. Eddy diffusion, here referring to turbulent processes that 
result from wave transience and dissipation, is parameterized in the model. To min- 
imize computational time required to calculate transport, all the short-lived species 
are grouped into families of species that interconvert rapidly. Then transport of the 
family a s a whole is calculated for the timestep. 

The boundary conditions for the chemistry part of the model include the 
ground-level sources and sinks of many gases such as N2O, CH4, CO2, and the 
chlorofluorocarbons. At each timestep, the chemistry calculations are carried out 
by making photochemical steady-state assumptions or by performing implicit or 
semi-implicit time integrations. Families are partitioned into member species. 

For a more detailed description of the model, see Crutzen and Gidel (1983). 


t 



CHAPTER 4 


USING TRACERS TO TRACK CHEMICAL 

CHANGE 

Our goal for this research is to follow chemical change in the stratosphere 
on timescales of weeks to years. In particular, we are looking at the response of 
relatively short-lived species — those with chemical lifetimes on the order of sec- 
onds to weeks — to perturbations such as an influx of sulfate aerosol. The species 
examined include ozone as well as members of two long-lived chemical families — 
the total inorganic chlorine family (Cly) and the reactive nitrogen family (NOy). To 
distinguish between chemical change and change due to transport in these species, 
we use the chemical tracers CH 4 and HF. Tracer analysis proves especially reliable 
when the species of interest has a short chemical lifetime but belongs to a family 
that is long-lived. As will be discussed below, in the case of ozone, tracer analysis 
is somewhat less trustworthy as a means to detect chemical change. 

Our work rests on the reasonable assumption that pairs of tracers show a 
linear correlation over a range of altitudes and latitudes (Kelly et al., 1989; Profitt 
et al., 1989; Profitt et al., 1990; and Plumb and Ko, 1992). Here tracers are defined 
as those species or families of species whose chemical lifetimes are very long relative 
to their lifetimes due to tranport. Lifetimes for transport processes or chemical 
loss are defined as the time it takes for a process to change the concentration of a 
species by a factor of 1/e. For chemical production, the lifetime is the ratio of the 
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concentration of the species divided by the rate of production. Two tracers that 
show linear correlation are said to be in gradient equilibrium, and their isopleths 
correspond to surfaces of rapid exchange of air parcels. The shapes of these rapid 
exchange surfaces are defined entirely by transport processes, independent of chem- 
ical sources and sinks. Gradient equilibrium between two species breaks down in 
regions where one of the species experiences a rate of chemical loss or production 
that is comparable in magnitude to the rate of transport processes. 

4.1. Background on the NOy and Cly Families 

and Ozone 

The two nitrogen species observed by HALOE, NO and NO 2 , belong to the 
NOy family. This family also includes HNO 3 , N 2 O 5 , CIONO 2 , HO 2 NO 2 , and a few 
other less abundant species. The source for nearly all NOy in the stratosphere is 
N 2 0, a gas emitted by biological and anthropogenic processes in the troposphere and 
carried to the stratosphere through the tropical tropopause. In the stratosphere, the 
main sink of N 2 0 is photolysis, which forms molecular nitrogen and atomic oxygen. 
But about 10 per cent of N 2 O reacts with 0( 1 D) to produce two NO radicals. The 
NOy family as a whole is long-lived throughout the stratosphere except at high 
altitudes (above about 45 km), 'where the photolysis of NO becomes an important 
NOy sink. The long chemical lifetime of the family means that transport processes 
play an important role in determining the distribution of NOy, and the abundance 
of NOy at 40 S is not invariant year after year, but may change with, for example, 
changing phase of the quasibiennial oscillation (QBO). We note that though the 
lifetime of the family as a whole is long, the lifetime of the member species against 



52 


conversion to another member of the family ranges from seconds (NO and N0 2 ) to 
weeks (HNO 3 ). Thus our research will focus on the chemical repartioning of the 
family in response to perturbations. 

The Cly family — which includes HC1, C10N0 2 , CIO, C10H, C1 2 0 2 , and Cl 
— is derived mainly from CH 3 CI, a natural tropospheric gas, and from anthropogenic 
sources — the chlorofluorocarbons (CFCs) and methyl chloroform (CH 3 CCI 3 ). As 
in the case of N 2 0, these source gases are transported upward through the tropi- 
cal tropopause into the stratosphere where they slowly photolyze to produce Cly. 
Transport back across the tropopause is the only sink of stratospheric Cly. Like 
NOy, Cly as a group has a long chemical lifetime in the stratosphere, on the order 
of years. Also like NOy, the member species have short lifetimes, ranging from 
seconds to weeks. 

Ozone in the stratosphere is produced via photolysis of molecular oxygen 
and destroyed via several catalytic cycles involving radicals such as NO or CIO 
(Wennberg et al., 1994). Production of ozone is greatest over the tropics where the 
flux of ultraviolet radiation is greatest. Loss processes dominate the ozone budget 
at high latitudes. 


4.2. THe Tracers CH 4 and HF 

Following the arguments of Plumb and Ko (1992), we may assume that 
parcels containing similar abundances of a tracer, such as CH 4 , also have similar 
abundances of the two long-lived chemical families Cly and NOy. This assumption 
allows us to follow the evolution of short-lived species that belong to the NOy and 
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Cly families — here NO, NOj, and HC1 — and to draw conclusions regarding the 
partitioning of Cly and NOy into member species even though not all the member 
species have been measured. 

HALOE measures the concentrations of two long-lived species — HF and 
CH 4 . Both species appear to be good candidates for our choice as tracer for this 
study. HF is derived entirely from chlorfluorocarbons (CFCs), which like the NOy 
source gas N 2 O are emitted in the troposphere and are slowly transported to the 
stratosphere. Production of HF is slow, with production lifetimes on the order of 
months to years, and transport back across the tropopause is the only sink. CH 4 , 
on the other hand, is itself emitted in the troposphere by both anthropogenic and 
natural processes. In the stratosphere, CH 4 slowly oxidizes by reaction with the 
radical species OH, 0( 1 D), and Cl. Both HF and CH 4 have chemical lifetimes on 
the order of years in the lower stratosphere and thus may be expected to correlate 
or anti- correlate linearly with other long-lived species or families. 

Figure 11 shows a plot of calculated NOy vs HF at various altitudes through 
the year at 40 S. The values were generated by the Max Planck time-dependent, two- 
dimensional, photochemical model (Bruhl and Crutzen, 1988; Prather and Rems- 
berg, 1993). Each loop represents 24 output points taken every 2 weeks throughout 

f 

the model year. The pressures represented range from 40 to 2 mb (about 20 to 40 
km). It is clear from the plot that up to about 800 ppt (10 mb), HF correlates well 
with NOy. Our task, therefore, will be to compare abundances of NOx in parcels 
that have similar abundances of tracer material — either HF or CH 4 . We may as- 
sume that such parcels contain similar abundances of NOy and that any change in 
NOx reflects changes in the partitioning of NOy and not changes in NOy itself due 
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to transport. Since NO and NO 2 have short chemical lifetimes against conversion 
to another NOy family member, we can assume that these species are sensitive to 
local conditions. 

Because the correlation between NOy and tracers breaks down in photo- 
chemically older air, we will exclude from our analyses air that contains more than 
800 ppt of HF. In other words, using HF as a measure of photochemical age, we 
can assume that signficant NO photolysis has taken place in air parcels of high HF 
content. In the mid-latitudes, air that has reached about 40 km and has spent about 
2 or 3 years in the stratosphere (Rosenlof, 1995) shows signs of the breakdown of 
gradient equilibrium between NOy and other tracers. 

A plot similar to Figure 11 may be constructed for NOy vs CH 4 , Cly vs HF, 
and Cly vs CH4. As noted above, Cly has no chemical sink in the stratosphere, and 
the correlation between Cly and tracers continues through high altitudes. Thus we 
can assume that parcels containing similar abundances of HF or CH 4 also contain 
similar abundances of Cly. Any changes in HC1 in these parcels over time must 
then be due to repartitioning of the Cly family. Like the NOy family members, 
Cly species have short lifetimes and can be expected to respond readily to local 
conditions. 

Tracking change in ozone poses special difficulties. In the regions focussed 
on in this study, ozone loss cycles dominate the ozone budget. These cycles pro- 
ceed more rapidly in air that is photochemically older and therefore contains greater 
abundances of radicals. Thus ozone anti-correlates with HF and correlates with CH 4 
in these regions. The long chemical lifetime of ozone, however, means that the re- 
sponse of ozone to perturbations is usually sluggish. While the NOy and Cly families 



55 


respond fairly rapidly to local conditions, ozone may take months. Nonetheless, as 
we will show below, tracer analysis can in some cases supply interesting information 
about changes in ozone concentrations due to chemical processes. 
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Calculated NOy vs HF through the year 
from 2 to 40 mb at 35 S 



HF (ppt) 


Figure 11. Annual cycles of NOy vs HF calculated by a two-dimensional, photo- 
chemical model for a range of altitudes (40 to 2 mb). Each loop represents a different 
altitude and connects the data points calculated every 15 days throughout the model 
year. 



CHAPTER 5 


RESPONSE OF SUMMERTIME ODD 
NITROGEN AND OZONE AT 17 MB TO 
MOUNT PINATUBO AEROSOL OVER THE 
SOUTHERN MID-LATITUDES 


Abstract 

Observations of stratsopheric NO, NO 2 , and O 3 from the Halogen Occulta- 
tion Experiment (HALOE) are examined over the southern mid-latitudes (35-45 S) 
in the summers following the eruption of Mount Pinatubo, from early 1992 through 
mid- 1995. The tracers HF and CH 4 , also observed by HALOE, are used to distin- 
guish between the effects of transport on the distribution of total reactive nitrogen 
(NOy) and ozone and the chemical effects of aerosol on NOy partitioning and ozone 
loss processes. HF mixing ratios are adjusted to take into account the increasing 
burden of stratospheric HF due to emissions of chlorofluorcarbons at the earth’s 
surface. Trends in aerosol surface area densities from the Stratospheric Aerosol 
and Gas Experiment II (SAGE II) are also studied for the same time period. Re- 
sults show that at 17 mb (about 28 km), in parcels of similar photochemical age, 
the abundances of summertime odd nitrogen (NO+NO 2 , or NOx) nearly doubled 
between 1992 and 1993 and then levelled off in subsequent years. The 1992-1993 
increase, coincident in time with a rapid drop in aerosol surface area density as 
the Pinatubo aerosol cleared the atmosphere, matches trends previously observed 
at lower altitudes (16-20 km). The trend demonstrates the importance of hetero- 
geneous chemistry in the partitioning of NOy even at these altitudes, where the 
higher photolysis rate of HNO 3 speeds up the cycling of nitrogen through the NOy 
family. Results also show that summertime ozone abundances, again in parcels of 
similar photochemical age, declined by nearly 10 per cent between 1992 and 1993 
and then, like NOx, remained about constant for the next three summers. The trend 
in ozone is opposite in sign to that observed at lower altitudes in the aftermath of 
the Pinatubo eruption and shows the importance of the NOx catalytic cycles in the 
ozone loss budget above about 25 km. 
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5.1. Introduction. 

The eruption of Mount Pinatubo in June 1991 injected about 20 million 
metric tons of S0 2 as high as 30 km into the stratosphere (Read et al., 1993; Bluth 
et al., 1992). On a timescale of weeks, the S0 2 was quickly converted to sulfuric acid, 
leading to surface area densities of sulfate aerosol 30 to 50 times that of background 
values (McCormick and Veiga, 1992; Deshler et al., 1993; Thomason et al., 1996). 
This large influx of aerosol perturbed the stratosphere in two important ways. First, 
the particles increased atmospheric heating by absorption of upwelling infrared ra- 
diation (Labitzke and McCormick, 1992), and in this way influenced stratospheric 
transport processes, including the timing of the phase change of the quasibiennial 
oscillation, or QBO (Kinne and Toon, 1992; Grant, 1995). Second, the particles pro- 
vided surfaces upon which heterogeneous reactions took place, leading eventually to 
changes in the abundances of key radicals in the ozone loss cycles. At mid-latitudes, 
the most important of these heterogeneous reactions is: 

N 2 0 5 (<7) + H 2 0(/) -» 2 HNO 3 G 7 ), (5.1) 

where the symbol g indicates that the species is in gas phase and /, liquid phase. 

Other heterogeneous reactions, such as the hydrolysis of C10N0 2 , achieve impor- 

< 

tance only at the cold temperatures present at high latitudes in winter (Hanson 
et al., 1994). 

Reaction 5.1 shifts nitrogen from N 2 Os, a relatively reactive form, to HNO 3 , 
a reservoir species with a lifetime on the order of weeks or months. Both N 2 Os 
and HNO 3 belong to the NOy family of molecules, which also includes NO, N0 2 , 
C10N0 2 , NO 3 , and H0 2 N0 2 . With the exception of HNO 3 , all the NOy molecules 
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are short-lived, and exchange of nitrogen between these molecules takes place on 
timescales of seconds to hours. An increased rate of Reaction 5.1 thus repartitions 
the NOy family, leading to decreased abundances of N0+N0 2 (NOx), which in 
turn, leads to changes in the abundances of other radicals and, in some regions, 
ozone. With less NO 2 available to sequester chlorine in C10N0 2> chlorine radical 
abundances increase, and chlorine-catalyzed ozone loss cycles speed up. Smaller 
abundances of NO, on the other hand, slow the reaction 

N0 + H0 2 -> N0 2 + 0H (5.2) 

and in this way shift the photochemical equilibrium between OH and HO 2 in favor 
of H0 2 . Because H0 2 is involved in the rate-limiting step of the odd hydrogen ozone 
loss cycle, more H0 2 results in a faster rate of ozone loss via that cycle. In the lower 
stratosphere, where the halogen and odd hydrogen loss cycles dominate the ozone 
loss budget (Wennberg et al., 1994), increased abundances of sulfate aerosol may 
thus indirectly lead to significant ozone losses. 

Above about 25 km, however, the NOx catalytic cycles which destroy ozone 
take on importance: 


I^O + O 3 — ► N0 2 + 0 2 
N0 2 + 0 — ► NO + 0 2 
net : O 3 + 0 — ► 2 0 2 . 


(5.3) 

(5.4) 


and 


NO + O 3 


N0 2 -f 0 2 


(5.5) 
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(5.6) 

(5.7) 


no 2 + o 3 -> no 3 + o 2 

N0 3 + hv - no + o 2 

net : 20 3 — ► 3 0 2 . 

On a per-radical basis, these cycles are less efficient than the halogen or odd hydrogen 
cycles, but with increasing altitude both NOx and the atomic oxgyen necessary for 
Reaction 5.4 become quickly more abundant. NOx therefore is predicted to be a 
dominant player in ozone removal processes at higher altitudes (Wennberg et al., 
1994). 

In the years following the eruption of Mount Pinatubo, stratospheric NOx 
loss coincident in time and space with the presence of Pinatubo aerosol was reported 
in many studies. Most of these studies, however, involved ER-2 observations made at 
low altitudes (Fahey et al., 1993; Kawaet al., 1993; Mills et al., 1993) or observations 
of total column amounts (Johnston et al., 1992; Coffey and Mankin, 1993; Koike 
et al., 1994). Webster et al. (1994) measured low abundances of N0 2 in the 20- 
30 km range near latitude 31 N on one day in August 1992. In the same height 
range, Rinsland et al. (1994) compared data from the March 1992 Atmospheric 
Trace Molecule Spectroscopy experiment to March 1979 data from the Limb Infrared 
Monitor and found elevated levels of HN0 3 in 1994. 

Ozone, on the other hand, showed losses in the lower stratosphere (12-22 
km) in 1992 and 1993, consistent with theory (Weaver et al., 1993; Hofmann et al., 
1994b). At altitudes above 24 km over North America during late 1992 and early 
1993, Hofmann et al. (1994b) observed higher than usual ozone concentrations and 
speculated that aerosol leftover from Pinatubo had suppressed NOx abundances 
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there, leading to enhanced ozone. Lacking simultaneous measurements of NOx or 
a tracer, however, Hofmann et al. could not rule out transport as a cause of the 
increased ozone. Measurements of total ozone from ground-based instruments or 
from the Total Ozone Mapping Spectometer (TOMS) on Nimbus 7 show significant 
decreases in column amounts in 1991 over the tropics and in 1992-1993 over the 
northern mid-latitudes (Gleason et al., 1993; Bojkov et al., 1993; Kerr et al., 1993; 
Herman and Larko, 1994; Randel et al., 1995). No significant decreases in total 
column amount were observed over the southern mid-latitudes. 

Data from HALOE provide for the first time a continuous record of NOx 
and ozone above 26 km in the aftermath of a major volcanic eruption. HALOE ob- 
servations of two tracers, HF and CH 4 , allow us to filter out the effects of transport 
on NOx abundances, including those transport effects brought on by the presence 
of Pinatubo aerosol itself. In the case of ozone, which has an intermediate chem- 
ical lifetime, transport effects are somewhat more difficult to filter out, as will be 
discussed below. We have focused this study on the southern mid-latitudes for two 
reasons. First, aerosol surface area densities in the Southern Hemisphere above 26 
km reached 30-40 times background values, while those in the Northern Hemisphere 
in the same altitude range peaked at only 8-9 times background values (Thomason 
et al., 1996; World Meteorological Organization, 1992). Second, aerosol loading in 
both hemispheres over the mid-latitudes peaked during late 1991 and early 1992 
— that is, during the summer months for the Southern Hemisphere and the winter 
months for the Northern Hemisphere. Given the shorter chemical lifetime of ozone 
during summer and the heavier aerosol burden of the Southern Hemisphere above 
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about 26 km, we expect to observe a more dramatic chemical response to Pinatubo 
aerosol in the summertime over the southern mid-latitudes. 

Probing the sensitivity of ozone to changes in NOx at these altitudes may 
provide some insight into the response of ozone to increased NOy emissions from 
the proposed fleet of high-speed civilian transports (HSCTs). Although the HSCTs 
are slated to fly at low altitudes in the stratosphere (about 18 to 20 km), it is not 
known to what extent the NOy emitted by the proposed fleet will be transported 
vertically. If, as has been suggested by Wofsy et al. (1994), vertical transport is more 
efficient than predicted by many current models, injection of NOy at 20 km may 
indeed affect NOx abundances at 28 km. While this paper examines the question 
of how decreases of NOx affect ozone abundances at high altitudes, the study may 
also shed light on how increases in NOx may affect ozone in the event that a fleet 
of HSCTs is launched. 


5.2. Method. 

Our goal for this work is to follow the evolution of NOx and ozone over 
the southern mid-latitudes at 17 mb (27-28 km) in the years following the Mount 
Pinatubo eruption — from summer 1992 through summer 1995. Figure 12 shows the 
surface area density of sulfate aerosol at 27.4 km at 40 S as a function of time from 
January 1991 through most of 1995 (Thomason et al., 1996). The righthand axis 
gives the surface area as a multiple of background surface area (World Meteorological 
Organization, 1992). The huge surge in surface area density in fall 1991 marks the 
arrival of the Pinatubo aerosol in this region. The trend is consistent with other 
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observations which show the bulk of the volcanic material confined to the tropical 
stratospheric reservoir in the first months following the eruption and then slowly 
dispersing to other latitudes (Trepte et al., 1993). At 40 S, the aerosol surface area 
density declined rapidly during 1992, relaxing back to pre-Pinatubo levels by the 
last months of the year. Though HALOE began observations in October 1991, the 
presence of huge abundances of aerosol made the early data unreliable at 17 mb. 
Useful measurements at this altitude were not made until January 1992. We must 
therefore look for evidence of the response of chemical species not to the injection 
of the Pinatubo aerosol, but to the gradual decline of the aerosol as the particles 
settled out of the atmosphere. 

It should be noted that lower in the stratosphere at this latitude, aerosol 
surface area densities remained at very high levels throughout 1992 and much of 1993 
(Thomason et al., 1996). Thus the data from lower altitudes, especially the NO and 
NO 2 measurements, show larger uncertainty due to the significant contribution of 
aerosol to extinction in the infrared. Even in aerosol-free conditions, however, the 
estimated total error in the measurements increases rapidly below about 25 km due 
to the smaller NO and NO 2 mixing ratios in that region. For example, the total 
error for NO 2 increases from about 10 per cent at 30 km to about 40 per cent at 
20 km (Gordley et al., 1996). fe ecause of the uncertainty in the data below 25 
km, especially under conditions of heavy aerosol, we do not analyze data from this 
region. 

To accomplish our task at 17 mb, we examine a series of similar air parcels 
through the years. In the case of NOx, this means we need to select parcels with 
similar NOy content so that changes in NOx will reflect changes in the partitioning 
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of NOy and not simply changes in the abundance of the NOy family as a whole. 
As described in Chapter 4, the NOy family as a whole throughout much of the 
stratosphere is long-lived and can be viewed as a tracer. Thus we can expect that 
the NOy family should correlate with the HALOE tracers HF and CH 4 . Changes in 
NOy partitioning can then be tracked by examining the NOx content of a series in 
air parcels that contain similar abundances of tracer material. Our calculations show 
that the timescale for conversion of NOx to HN0 3 is about one week, consistent with 
the observations of Kawa et al. (1993). We may infer then that NOy partitioning 
should indeed be sensitive to the local aerosol loading. 

As noted in Chapter 4, ozone anti-correlates with HF in this region due to 
the greater abundance of radicals in photochemically older air. Unlike in the case 
of NOx, however, 0 3 at this altitude has a chemical lifetime of about 70 days in 
summer, comparable to its lifetime due to transport (Brasseur and Solomon, 1986; 
Ko et al., 1989). Since the ozone at any given point in space will respond slowly to 
the local conditions while retaining a “memory” of conditions at another point, we 
expect therefore only weak anti-correlation or correlation of ozone with the tracer 
gases HF or CH 4 . Nonetheless, as will be shown below, using tracers in our analysis 
will allow us to define with greater certainty the ozone response to the Pinatubo 
aerosol. 

The next question to address in this study is which of the two tracers 
measured by HALOE best serves our purposes. Figure 13 shows a plot of observed 
CH 4 vs HF at 17 mb over the latitude range 35 to 45 S duriug the time interval 
from early 1992 through most of 1995. In this part of the atmosphere, the lifetimes 
of both HF and CH 4 are very long, on the order of years, so one would expect the 
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two species to exhibit a compact linear relationship (Plumb and Ko, 1992). Since 
the stratospheric burden of HF is increasing rapidly due to increasing emissions of 
CFCs at the earth’s surface, we have adjusted the HF mixing ratios in Figure 13 as 
follows. The observed value of HF on any given day was multiplied by the factor 
(1 + 0.066i /36-5) — 1 where t is the time in days since the start of the observations and 
6.6%/year is the annual growth rate of stratospheric HF (Russell et al., 1996). We 
have not adjusted the CH 4 mixing ratios since the annual increase for global CH 4 
is less than one per cent (Dlugokencky et al., 1994). The correlation coefficient for 
CH 4 vs the adjusted HF is 0.91 for the 1280 data points, compared to only 0.81 for 
CH 4 vs unadjusted HF. The improved correlation coefficient for the adjusted data 
gives us confidence in our correction. 

All of the scatter that remains in Figure 13 after the HF adjustment can 
be assigned to random instrument and retrieval algorithm error, estimated to be 
about 7 per cent for both CH 4 and HF at this altitude. We note that 7 per cent 
may, in fact, be an overestimate of the random error for these channels, given the 
extremely good correlation of the two species. In any case, since our work rests on 
the assumption that NOy as a family correlates well with long-lived tracers, we have 
chosen for this study to look at only those parcels whose HF-CH 4 relationships lie 
with two standard deviations of the least squares fit of the data in the Figure, about 
55 per cent of the points. In this way, we use information from both tracers in our 
analysis. 

We note here that for the purposes of this study, we have assumed that 
the eruption of Mount Pinatubo did not inject appreciable quantities of fluorine into 
the stratosphere. Although Symonds et ah (1988) predict that explosive volcanic 
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eruptions could transmit HC1 and HF to the stratosphere, the model of Tabazedeh 
and Turco (1993) shows that HC1 vapor in a volcanic plume is efficiently scavenged 
by the huge quanities of water droplets also present in the plume. The study of 
Tabazedeh and Turco did not include HF. HF is somewhat less soluble in water 
than is HC1, but we detect no perturbations in HF-CH 4 correlation during the time 
period of our study. Thus, we have assumed that, like HC1, whatever HF was 
emitted during the eruption was scavenged by water droplets, followed by rainout. 

5.3. Observations 

Figure 14 shows the relationship of measured summertime NOx to HF at 
17 mb over the latitude region 35 S to 45 S for the four years following the eruption 
of Mount Pinatubo. For this and subsequent plots, the HF data have been adjusted 
to take into account the increasing emissions of CFCs, and all data points with 
CH 4 -HF relationships outside two sigma of the least squares fit of Figure 13 have 
been filtered out. The red triangles represent data from 1992, the dark blue from 
1993, the light blue from 1994, and the green from 1995. The 1993 data include 
some points from the last few days of December 1992. It is clear from Figure 14 
that for similar values of HF, summertime NOx was markedly lower in 1992 than in 
subsequent years. 

Note that plotting NOx vs HF allows us to distinguish easily between 
parcels with low NOx due to low NOy and parcels with low NOx due to unusual 
partitioning of NOy. Parcels with low HF values have spent less time in the strato- 
sphere and have been exposed to less ultraviolet radiation than parcels with high 
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HF. Such parcels almost certainly have low NOy abundances, and consequently low 
NOx. Low NOx values in parcels with high HF, on the other hand, signify unusual 
partitioning of NOy towards more HNO 3 . It is worth noting that in summer 1994, 
no parcel with an HF value less than about 480 ppt was observed, while during the 
other three summers, parcels with as little as 350 ppt were detected. This implies 
that in 1994 NOy mixing ratios were also higher than usual for reasons most likely 
having to do with transport. 

In Figure 15a, we show the evolution of summertime NOx through the 
years for parcels characterized by an HF content of 500 to 550 ppt — that is, for 
parcels whose NOy content also falls within a similarly narrow range. Each data 
point is represented by a dot, while the average NOx and date for each summer is 
denoted by a triangle. The solid line connects the summer averages, and the dotted 
lines show the range of two standard deviations from the mean of the last three 
summers, 1993-1995. Not surprisingly, the plot shows extremely low NOx values 
for 1992, followed by a levelling off in subsequent years to almost double the 1992 
values. The scatter in Figure 15a may be attributed to several sources. First, the 
HF increment of 50 ppt corresponds to an increment of about 0.3 ppb of NOx, so 
some scatter in the plot is to be expected. Second, the partitioning of NOy may 

‘r 

vary somewhat with latitude due to differences in solar insolation, leading to a range 
of values of NOx for the same HF value. Lastly, instrument error almost certainly 
contributed to the scatter. Despite the scatter, however, the trend in NOx is clear, 
and all the points from 1992 lie outside the two-sigma range of the mean of the three 
subsequent years. 
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Figure 16 shows summertime 0 3 vs HF for the same time interval and 
latitude range as Figure 14. The HF data have been adjusted as in Figure 14, and 
the colors of the triangles correspond to the same set of years. The correlation 
coefficient R for the 125 data points in the plot is -0.5, which indicates that the 
two species are only weakly anti-correlated, as expected. As noted above, the anti- 
correlation comes about since photochemically older air contains greater abundances 
of the radical species that destroy ozone. Much of the scatter in the plot probably 
arises from the long chemical lifetime of 0 3) which means that the gas may retain 
evidence of chemical processes occurring in other regions. Nonetheless, we see that 
ozone abundances from the 1992 summer are among the highest for any given value 
of HF. 

As in Figure 15a, we plot in Figure 15b the summertime 0 3 values for 
those parcels whose HF content ranges between 500 to 550 ppt. As before, the dots 
represent all the data points, the triangles the season averages, and the two dotted 
lines signify the range of two standard deviations from the mean mixing ratios for 
the 1993-1995 summers. The mean 0 3 mixing ratio for summer 1992, about 10 per 
cent greater than the mean of the following three summers, lies just outside the 
two-sigma range, indicating that ozone abundances were indeed unusual that year. 

As a check on our methods, we have replotted Figures 15a and 15b, making 
use of all the summertime data available. In Figures 17a and 17b we have not 
imposed any limits on HF content of the points included; all data points for each 
summer at 17 mb are shown. Considerably more scatter is apparent in both plots. 
In Figure 17a, we see that NOx again increases between summer 1992 and summer 
1993, but then continues to increase at nearly the same rate during the subsequent 
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year. Summertime NOx in 1995 shows a small drop. We believe that the high NOx 
data of 1994 merely reflect the high NOy values of that year, as predicted by the 
high HF values of Figure 14. In other words, the 1994 change in NOx is probably 
due to transport processes and not to chemistry. In Figure 17b, the dotted lines 
again show the two-sigma range for the 1993-1995 data, but this time summertime 
ozone for 1992 falls inside the range. The slight increase in ozone in 1995 may be 
linked to the decrease in NOx for that year. Since our goal is to determine changes 
in NOx and O3 due solely to chemistry, Figures 17a and 17b are not useful to our 
study. As will be discussed in the next section, the trends in the unfiltered data are 
not consistent with the observed trends in surface area densities of sulfate aerosol. 

A pair of plots similar to Figures 15a and 15b may be obtained for HALOE 
data from the latitude region 25 S to 35 S, again at 17 mb. In this region of 
photochemically younger air (compared to 35-45 S), the HF interval chosen was 450 
to 500 ppt. In these plots (not shown) summertime NOx again nearly doubles from 

1992 to 1993, and O3 decreases by about 7 per cent. The smaller ozone change for 
the 25-35 S latitude range, compared to that for the 35-45 S latitude range, may 
reflect the lesser importance of ozone loss processes over lower latitudes. 

5.4. Discussion and Model Calculations. 

The dramatic rise in NOx abundances between summer 1992 and summer 

1993 coincident in time with the rapid decline in aerosol surface area density demon- 
strates the importance of Reaction 5.1 in the partitioning of the NOy family at 17 
mb in the mid-latitude stratosphere. The data suggest that by early 1992, when 
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HALOE observations began in this region, the presence of large abundances of sul- 
fate aerosol leftover from the Mount Pinatubo eruption had led to a shift in reactive 
nitrogen from NOx to the reservoir molecule HN0 3 , which even at these altitudes 
did not photolyze quickly enough to maintain NOx at its usual level. As the aerosol 
cleared the atmosphere, NOx concentrations apparently returned to pre-Pinatubo 
values. 

We note that HNO3 data from the Cryogenic Limb Array Etalon Spec- 
trometer (CLAES), another instrument aboard UARS (Roche et al., 1993), agree 
well with this interpretation of the trends in HALOE odd nitrogen. CLAES ceased 
observations in April 1993, so only two summers of data are available for a compar- 
ison with HALOE. Zonally averaging CLAES HN0 3 at the 17 mb level between 42 
and 44 S for the same days as the HALOE observations yields an average 6.6 ppb 
for summer 1992 and 4.1 ppb for summer 1993. This drop in mixing ratio of 2.5 ppb 
is larger than the corresponding increase in NOx by about 1.2 ppb. The missing 
1993 NOy was probably tied up in N2O5, which during the night builds up to higher 
values when more NOx is available at sunset. Note that the HALOE observations 
are taken at sunrise, when N 2 Os mixing ratios reach their peak. 

The response of ozone to the increased burden of sulfate aerosol at this 
altitude was more complicated. As noted in the Introduction, NOx plays both a 
direct and an indirect role in the ozone budget. The direct role involves the two 
NOx catalytic cycles that destroy ozone, Reactions 5. 3-5.7. The indirect role involves 
the sequestering of active chlorine by N0 2 in the molecule C10N0 2 as well as the 
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competition between the Reaction 5.2 and the rate-limiting step (Reaction 5.9) in 
the odd hydrogen catalytic cycle: 

0 3 + 0H -► O 2 + HO 2 (5.8) 

O 3 4- HO 2 — > 2 O 2 + OH (5.9) 

fict : 2 O3 — ♦ 3O2. 

NOx thus controls the abundance of chlorine radicals and the rate of Reaction 5.9. 
While the direct role of NOx in the ozone budget leads to the removal of ozone, the 
indirect role leads to slower rates of destruction of ozone by active chlorine and odd 
hydrogen radicals. At altitudes probed by the ER-2 (16-20 km), the halogen and 
odd hydrogen cycles dominate ozone removal processes, while the NOx catalytic 
cycles account for only about 20 per cent of ozone loss (Wennberg et al., 1994). 

The significant loss of NOx in the bottom few kilometers of the stratosphere that 

followed the eruption of Mount Pinatubo thus led to depressed ozone concentrations 
in that region (Weaver et al., 1993). 

A few kilometers higher than the ER-2 observations, however, the response 
of ozone to perturbations in NOx abundances changes. Although the NOx catalytic 
cycles are less efficient than the other ozone loss cycles, NOx abundances increase 
rapidly with height because HNO 3 photolysis rates increase with the increased solar 
insolation available at higher altitudes. NO 3 photolysis rates also increase at higher 
altitudes, speeding up Reaction 5.7, and atomic oxygen becomes more abundant, 
resulting in a faster rate of Reaction 5.4. Both Reaction 5.4 and Reaction 5.7 are 
the rate-determining steps in their cycles. Although an important direct role for 
NOx in the 0 3 budget has been predicted at the higher altitudes (e.g., Wennberg 
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et al., 1994), no unambiguous evidence of the sensitivity of 0 3 to NOx above about 
25 km has previously been observed. Figure 15b shows for the first time the sign of 
the change in ozone at 17 mb that accompanies a loss in NOx at that altitude. The 
plot suggests that although the suppressed levels of NOx in summer 1992 may have 
increased the rates of 0 3 removal via the active chlorine and odd hydrogen cycles, 
they also decreased the rates of 0 3 loss via the NOx catalytic cycle, leading to a net 
enhancement in 0 3 mixing ratios for this altitude. 

To test these ideas, we have performed calculations with a time-dependent, 
photochemical box model described in Chapter 3. We initialized the model with 
data from HALOE and ran a series of simulations at two altitudes (28 and 20 km), 
varying the aerosol surface area densities from background conditions to 50 times 
background conditions. For each run, initial mixing ratios were input at the sunrise 
timestep, and the model was run to diel steady-state, as described in Chapter 3. 
Data for NO, NO2, 0 3 , HC1, H2O, and CH4 were taken from HALOE observations 
on February 1 1 , 1992, when the mean latitude of HALOE observations was 43 S. For 
observations at 28 km (17 mb), we took the average of those data points whose HF 
content fell within the range of 500 to 550 ppt. For observations at 20 km (54 mb), we 
took the average of all the data points on that day, keeping in mind that the aerosol 
had a larger effect on NOx measurements in this region than at higher altitudes. 
We determined values for Cly and NOy at both altitudes by first calculating N 2 0 
from HALOE CH4, making use of a relationship between these two tracers derived 
from aircraft data (Woodbridge et al., 1995). The abundances of Cly and NOy were 
then derived from that of N 2 0 by the methods of Woodbridge et al. (1995) and 
Fahey et al. (1990), respectively. An initial value of C10N0 2 was determined by 
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simply subtracting HALOE HC1 from Cly. We assumed N 2 0s to be negligible at the 
high surface area densities of February 1992, and after C10N0 2 and NOx had been 
accounted for, we assigned the remaining NOy to HN0 3 . In fact, since the model was 
run to diel steady-state conditions, the output was nearly insensitive to the input 
concentrations of the individual Cly and NOy member species (though not to the 
input values for the total families). For the bromine species, we used 5 ppt for BrO 
and 10 ppt for total inorganic bromine, or Bry (Avallone et al., 1995). Background 
aerosol densities were taken from the World Meteorlogical Organization (1992). To 
calculate the appropriate radiative flux and photolysis rates at the two altitudes 
studied, ozone column amounts were taken from HALOE data. A summary of key 
input values for the two altitudes is shown in Table 1. 

As an aside, we note that we chose not to initialize the model with data 
from CLAES. Observations from CLAES and HALOE were not coincident in time 
and space. Zonally averaging CLAES N 2 0 data from February 11, 1992, for the 
latitude band 42-44 S yields a value of N 2 0 not consistent, according to the “rule” 
of Woodbridge et al. (1995), with the CH4 mixing ratios observed by HALOE. 

Figure 18 shows calculated NOx as a function of the multiple of background 
aerosol surface area density. The solid line represents data from 28 km, and the dot- 
ted line from 20 km. At each altitude, the concentrations have been normalized to 
the abundance of NOx calculated for background conditions of sulfate aerosol. As 
the surface area density of sulfate aerosol increases from 1 to 10 times background 
conditions, NOx abundances at both altitudes decline sharply, signifying a signifi- 
cant shift from NOx to HNO3 in the presence of high aerosol concentrations. The 
slightly less rapid response of NOx at 28 km to increasing sulfate aerosol is caused 
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by the higher photolysis rate of HN0 3 at that altitude (about double the rate at 
20 km). As expected, HNO3 is a slightly less stable reservoir for NOy at the higher 
altitude. 

At surface area densities greater than about 10 times background, NOx 
abundances begin to level off, a sign of the saturation of N2O5 hydrolysis (Reac- 
tion 5.1) on sulfate aerosol (Mills et al., 1993). Saturation is achieved when the 
photolysis rate of N 2 0 5 becomes much slower than the rate of N 2 0 5 hydrolysis, and 
the rate-limiting step in the heterogeneous conversion of NOx to HNO3 becomes 
the three-body reaction of N0 2 and N0 3 to form N 2 0 5 , a process not dependent 
on surface area of aerosol. In other words, under conditions of high aerosol, the 
photolysis rate of N 2 0 5 cannot compete with the rate of hydrolysis, and during the 
day little nitrogen is cycled back to NOx. This means that insufficient N 2 Os is avail- 
able to keep the heterogeneous reaction going as rapidly as might be expected. At 
night, when photolysis of N 2 0 5 ceases, the concentration of that molecule builds. 
With increasing sulfate aerosol, however, less NOx becomes available at night to 
form N 2 0 5 , and the nighttime hydrolysis of N 2 Os on aerosol proceeds more and 
more slowly. As shown in Figure 18, the saturation of N 2 Os hydrolysis is achieved 
at a lower surface area density at 20 km than at 28 km because of the slower N 2 0 5 
photolysis rates at the lower altitude. 

According to Figure 12, SAGE II measured surface area densities of sulfate 
aerosol at about 25 times the background value in this region of the atmosphere in 
February 1992. The calculated per cent change in NOx that accompanies a change 
in sulfate aerosol from background conditions to 25 times background conditions is 
about 60 per cent, consistent with the observations in Figure 15a. 
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In Figure 19, as in Figure 18, we plot O3 as a function of multiple of aerosol 
surface area density. Again we show data from the two altitudes and normalize the 
data at both altitudes to the concentration of ozone under background conditions 
of sulfate aerosol. Calculated ozone abundances at 20 km decline with increasing 
aerosol surface area, consistent with earlier observations. At 28 km, calculated ozone 
concentrations increase with increasing aerosol surface area, reaching a value about 8 
per cent greater than the value under background conditions. Above about 25 times 
background levels of aerosol, ozone concentrations begin to level off. The calculated 
ozone increase at 25 times background is consistent with the model results of Tie et 
al. (1994) in this region and is slightly less than the observed ozone increase of 10 
per cent. 

To gain understanding of the behavior of O3 at 28 km, we plot in Figure 20 
the rates of the major ozone loss cycles, grouped by family, again as a function of 
multiple of surface area density. The dotted line represents the total loss rate due to 
the halogens, the dashed line that due to odd hydrogen, and the solid line that due to 
NOx. The loss rates, in units of molecules / cm? / s, have been normalized to the loss 
rate of NOx at background surface area densitites. For each family, the rates of the 
rate-limiting steps of each important cycle have been summed together. All rates are 
for February at 28 km at 40 S. At low multiples of background aerosol, NOx clearly 
dominates the ozone loss cycles, with rates about three times that of the halogen 
cycles. With increasing aerosol concentrations, NOx loss rates first decline and 
then level off, consistent with the levelling off of NOx abundances under saturation 
conditions. At higher aerosol surface areas, with less NOx available to sequester 
CIO in C10N0 2 , the halogen cycles gain in importance in the ozone budget. The 
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faster loss rates of the halogen loss cycles, together with the levelling off of the NOx 
loss rates, limit to some extent the ozone increase that accompanies an increase in 
the sulfate aerosol burden at this altitude. Thus in Figure 19, as aerosol surface area 
density increases from 30 to 50 times background, the calculated ozone abundance 
at 28 km changes remarkably little. 

The different responses of ozone to perturbations in NOx at the two alti- 
tudes studied give insight into the trends observed in total ozone in the years fol- 
lowing the eruption of Mount Pinatubo. In the Southern Hemisphere, where aerosol 
loading at altitudes above about 25 km reached 30-40 times background (Thomason 
et al., 1996; World Meteorological Organization, 1992), significant ozone enhance- 
ment took place and may have offset the decreases in ozone at lower altitudes. This 
may explain why total column amounts of ozone in the Southern Hemisphere showed 
little or no change in the aftermath of the Pinatubo eruption. In the Northern Hemi- 
sphere, smaller aerosol abundances perturbed the region above 26 km, and the total 
column amounts of ozone reflect mainly the losses at lower altitudes. 

5.5. Summary and Implications. 

HALOE observations in the southern mid-latitudes of summertime NOx 

f 

and ozone, together with SAGE II measurements of aerosol surface area densities, 
show the sensitivity at 17 mb of NOy partitioning and ozone loss processes to the 
presence of huge abundances of aerosol leftover from the 1991 eruption of Mount 
Pinatubo. As the aerosol cleared this region of the atmosphere, summertime NOx 
concentrations nearly doubled between 1992 and 1993, and then held steady through 
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summer 1995. The trend in NOx matches trends observed at lower altitudes and 
shows that even at 17 mb, where HN0 3 photolyzes twice as rapidly as in the lower 
stratosphere, NOy partitioning is extremely sensitive to the presence of aerosol. 
Ozone, on the other hand, decreased by about 10 per cent from summer 1992 to 
summer 1993, suggesting that the NOx catalytic cycles play a major role in the 
ozone budget at these altitudes. In the aftermath of the Pinatubo eruption, the loss 
of NOx, while leading to accelerated rates of the halogen and odd hydrogen loss 
cycles in the ozone budget, also resulted in slower rates of the NOx loss cycles. The 
net effect at this altitude, where the NOx cycles dominate the ozone budget, was a 
temporary increase in ozone. As NOx returned to pre-Pinatubo values, the rates of 

the NOx catalytic cycles sped up, and ozone abundances decreased and then levelled 
off. 

These trends in NOx and 0 3 were determined by using the tracers HF and 
CH 4 , two other molecules observed by HALOE. Use of these tracers enabled us 
to compare parcels of similar photochemical age, containing similar abundances of 
NOy. For NOx, this meant that any changes we observed in N0+N0 2 abundances 
were due entirely to changes in the partitioning of NOy, and not due to trans- 
port effects on NOy itself. For 0 3 , this meant we could compare air parcels which 
under normal conditions would contain similiar abundances of radicals, leading to 
similiar rates of ozone loss. Ozone at this altitude anti-correlates with HF due to 
the increasing abundance of radicals in photochemically older air. Although this 
anti-correlation is somewhat weak due to the long chemical lifetime of ozone, tracer 
analysis allowed us to discern a trend in ozone which is in good agreement with the 
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trend predicted by our photochemical bo* model. Thus tracer analysis has proven 
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Table 1. Initial conditions for photochemical box model 
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Monthly averaged aerosol surface area density 
at 27.5 km and 40 S 
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Sunrise CH 4 vs HF at 17 mb 
from late 1991 - late 1995, 35 S to 45 S 
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Figure 13. HALOE CH 4 vs HF from early 1992 through most of 1995 on the 17 
mb surface over the southern mid-latitudes (35-45 S). HF has been adjusted to take 
into account the annual increase of stratospheric HF due to increasing emissions of 
CFCs at the earth’s surface (see text). 
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Summertime NOx vs HF at 17 mb 
between 35 and 45 S, 1992-1995 



Figure 14. Summertime NO+NO 2 vs IIF at, 17 mb between 35 and 45 S during 1992- 
1995. The color for each year is indicated in the legend. HF has been adjusted, and 
the data filtered, a,s described in the text. 
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Figure 15. Summertime (a) NO+NO 2 and (b) O3 at 17 mb over the southern 
mid-latitudes during the years 1992 through 1995. Only those points representing 
parcels with an HF content between 500 and 550 ppt are shown. Each dots denotes 
a separate observation, and the triangles show the seasonal mean for that year. The 
solid line connects the seasonal means, and the two dotted lines show the two-sigma 
range of the mean of the 1993-1995 summers. 
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Summertime 03 vs HF at 17 mb 
between 35 and 45 S, 1992-1995 
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imertime O 3 vs HF at 17 mb between .35 and 45 S dur 
4, the color for each year is indicated in the legend, 
he data filtered. 
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NOx and Ozone from 1992 to 1995 
summer 35 S to 45 S 



Figure 17. As in Figure 15. Here all data from parcels on the 17 mb surface, 
regardless of HF content, have been plotted. 
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Normalized NO+NO2 at two Altitudes 



Figure 18. Calculated NOx vs multiple of surface area density. The solid line 
represents model output at 28 km; the dotted line at 20 km. 
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Normalized O3 at two Altitudes 



Figure 19. Same as Figure 18, but for O 3 . 
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Rates of Ozone Loss Cycles 
vs Multiple of Background Surface Area 



Multiple of background surface area 


Figure 20. Calculated rates of the major ozone loss cycles, grouped according to 
family, vs multiple of background surface area density at 28 km. For each family, 
the rates of the rate-limiting steps for each cycle have been summed together. The 
solid line denotes the NOx cycles, the dotted line halogen cycles, and the dashed 
lines the odd hydrogen cycles. The loss rates, in units of molecules/ cm 3 / s, have 
been normalized to the loss rate of NOx at background surface area densitites. 




CHAPTER 6 


EVOLUTION OF CHLORINE AND NITROGEN 
SPECIES IN THE LOWER STRATOSPHERE 
DURING ANTARCTIC SPRING 


Abstract 

Observations of O 3 , HC1, NO, and N0 2 from the Halogen Occultation Ex- 
periment (HALOE) provide a means to investigate chemical change in the lower 
stratosphere over Antarctica during 23 days in October of 1992. Two long-lived 
species also observed by HALOE, HF and CH 4 , are used as tracers to identify 
a series of air parcels having similar Cly and NOy abundances and similar tem- 
perature histories. Most of the parcels chosen using tracer analysis show low O 3 
mixing ratios, less than 1 ppm on the 480 K surface (about 20 km). Parcels with 
O 3 mixing ratios higher than 1 ppm were occasionally observed in regions of steep 
horizontal gradients of tracer concentrations where mixing with mid-latitude air 
may have occurred. HC 1 mixing ratios for those parcels with less than 1 ppm 
of ozone more than doubled from early to late October, and NO and N0 2 abun- 
dances rose sharply, by a factor of five for NO and a factor of three to four for 
N0 2 . These trends in HC1, NO, and N0 2 agree with model calculations which 
show that the formation of HC1 proceeds quickly when O 3 levels fall so low that 
( 1 ) the rate of the reaction Cl + O 3 — * CIO + 0 2 slows and ( 2 ) the rate of the 
reaction NO + CIO — * N0 2 + Cl becomes faster than the rate of the competing re- 
action NO -f O 3 — ♦ N0 2 + 0 2 . Under these conditions, Cl increases at the expense 
of CIO, and HC1 is formed via the reaction Cl + CH 4 — ► HC1 + CH 3 . Stratospheric 
chlorine is thus shifted from re'hctive species to the long-lived, reservoir molecule 
HC1. The repartitioning of the active chlorine family in favor of HC1 may limit the 
extent of ozone depletion. In addition, as the chlorine in C10N0 2 is shifted to HC1 
via CIO and Cl, more active nitrogen becomes available to form NO and N0 2 . 
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6.1. Introduction 

In recent years, in-situ and satellite measurements, together with labo- 
ratory experiments and model calculations, have increased our understanding of 
the chemical and transport processes leading to the rapid depletion of stratospheric 
ozone each spring over Antarctica. (See, for example, Solomon et al., 1986; Solomon 
et al., 1990; and the special issues of the Journal of Geophysical Research, 94, nos. 
D9 and D14). The conditions for the annual springtime ozone losses are set the 
previous winter, when the steep temperature gradient between the pole and the 
equator leads to strong circumpolar winds, or the polar night jet. The winds act as 
a barrier to limit horizontal mixing between polar and mid-latitude air. As temper- 
atures fall during winter, polar stratospheric clouds (PSCs) form, providing surfaces 
upon which heterogeneous chemical reactions may take place. These reactions shift 
reservoir forms of stratospheric chlorine — HC1 and C10N0 2 — to photochemically 
active forms, such as Cl 2 or C10H. The reactions also shift stratospheric nitrogen 
from active forms to the reservoir species HNO 3 . In spring, when increased sunlight 
returns to the region, the active chlorine species readily photolyze, leading to large 
abundances of Cl, which then forms CIO via Cl + O 3 — ► CIO + 0 2 . The radicals 
Cl and CIO both participate in several catalytic cycles that rapidly destroy ozone. 
Since active nitrogen abundances are low, there is little N0 2 available to react with 
CIO to form the less reactive species C10N0 2 , and the chlorine catalytic cycles 
continue unimpeded. 

While much recent work has studied the processes leading to ozone loss, 
less attention has been paid to the transport processes or chemical mechanisms that 



91 


shut off the fast catalytic cycles that destroy ozone. Mixing of perturbed with unper- 
turbed air not only disperses the chlorine radicals, but also permits these radicals to 
react with the larger abundances of active nitrogen species present in unperturbed 
air. These reactions transfer chlorine to longer-lived, less reactive species. Schoeberl 
et al. (1992) have used aircraft data and model calculations to show that mixing is 
slow within the vortex, and that mixing at the vortex edge is essentially one-way, 
with filaments of vortex air being stripped away toward mid-latitudes. More re- 
cently, Schoeberl et al. (1995) have analyzed HALOE data and found evidence that 
the vortex itself remains isolated from mid-latitude air through October. There is, 
however, much evidence that suggests that rapid mixing between midlatitude and 
vortex air occurs in a band about 10 degrees latitude in width, centered near the 
polar night jet (e.g., Pierce et al., 1994a, and Tuck et al., 1995). 

Regardless of the extent of mixing with midlatitude air, chemical processes 
within the vortex may put a brake on the chlorine catalytic cycles once ozone abun- 
dances drop to a certain threshold mixing ratio. The theoretical work of Prather 
and Jaffe (1990) predicted that in typical PSC-processed air — that is, air char- 
acterized by low abundances of ozone and high abundances of active chlorine — 
chlorine may suddenly shift from active forms to the long-lived species HC1 through 

t 

reactions involving NO and CH4. This quick shift of chlorine to a reservoir species 
effectively shuts down the ozone loss cycles. 

Using data from UARS, Douglass et al. (1995) have reported evidence of 
this shift in the Antarctic, where springtime ozone levels drop to very low levels, 
but not in the Arctic, where ozone is less depleted. Douglass et al. (1995) used 
modified potential vorticity (MPV) to bin parcels for their analysis. This approach 
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is based on the assumption that parcel trajectories tend to follow potential vorticity 
contours, and that parcels with the same potential vorticity have experienced similar 
conditions, such as time exposed to polar stratospheric clouds. Some recent studies, 
however, suggest that in the lower stratosphere inside the Antarctic and Arctic 
vortices, potential vorticity may not correlate as well with long-lived tracers such as 
N 2 O as has been assumed (e.g., Tao and Tuck, 1994; Yang, 1995). Diabatic descent 
or mixing at the vortex edge may account for the deviations from the expected linear 
correlation between potential vorticity and long-lived tracers over time. As will be 
discussed below, these deviations pose difficulties in monitoring the repartitioning 
of long-lived chemical families into member species. 

In our work, we have used the long-lived tracers CH 4 and HF to identify 
parcels within the Antarctic vortex that are of similar photochemical age, have 
experienced similar rates of denitrification, and thus have similar abundances of 
both total inorganic chlorine (Cly) and active nitrogen (NOy). In the series of air 
parcels chosen by tracer analysis, we detect a rapid rise in HC1 during the spring of 
1992, accompanined by similarly rapid rises in NO and NO 2 . These trends are in 
qualitative agreement with those detected by Douglass et al. (1995). We have also 
used HALOE tracers CH 4 and HF, together with data from the Cryogenic Limb 
Array Etalon Spectrometer (CLAES), to determine the Cly and HNO 3 content of 
the parcels under investigation. Using this information and a photochemical box 
model, we have modeled the time evolution of HC1, NO, and NO 2 . 
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6.2. Use of Tracers to Select Air Parcels 

In 1992, HALOE sunrise observations swept toward and then away from 
the southern polar region in the last half of September. HALOE sunset observations 
then moved close to the polar region and then away during the month of October. 
We have used this sequence of measurements in our study of the chemical evolution 
of trace species during Antarctic spring. 

In order to follow the evolution of short-lived chemical species in this vor- 
tex air using HALOE data, it is essential to identify a series of parcels of similar 
photochemical age and with similar histories — that is, parcels that contain like 
amounts of total inorganic chlorine and active nitrogen and that have experienced 
comparable temperatures. For this work, we have used long-lived chemical tracers 
to select a series of similar parcels during a period of 23 days beginning October 1 , 
1992. 

Chapter 4 outlines the strategy for using tracers such as HF and CH 4 
for our work. Briefly, we showed there that for a series of air parcels containing 
similar abundances of the tracers HF or CH 4 , we could assume that the parcels also 
contained similar abundances of the families NOy and Cly. In this way, we could 
assign any changes observed in .the member species of these families over time to 
chemical processes that repartition the family. 

We note here that within the polar vortex, parcels of air may experience 
NOy loss through the process of denitrification, in which cloud particles containing 
HNO 3 grow to sufficient size to fall to lower altitudes. For our work, we assume 
that mixing takes place rapidly along isopleths, and thus parcels containing similar 
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quantities of a tracer such as CH 4 have also experienced similar temperature histories 
and similar growth rates of polar stratospheric clouds. These conditions would lead 
to similar rates of denitrification for all the parcels under consideration, and similar 
abundances of NOy leftover after denitrification. 

Much recent work has documented the descent of photochemically older 
air within the Antarctic vortex during the winter months (e.g., Russell III et al., 
1993b). Such air is characterized by low abundances of CH 4 and high abundances 
of HF. Figure 21 shows the distribution of sunrise HF on the 480-K surface during 
September and early October of 1992. High values of HF mixing ratio (greater than 
500 ppt) are apparent in the region poleward of 60 S in the region between about 
90 W and 90 E. The mid-latitudes northward of about 45 S are characterized by 
mixing ratios of 200 ppt or less. In between these two extremes, a steep horizontal 
gradient in mixing ratios is observed. 

Figure 22 shows the correlation of HF with CH 4 at high latitudes in the 
lower stratosphere during the period October 1-23, 1992. A large degree of scatter 
is apparent in this plot. There are three possible sources for the scatter in Figure 
22. First, the annual increase in stratospheric fluorine due to increasing emissions at 
the earth’s surface may affect the correlation between HF and another tracer when 
air characterized by a range of ages is plotted. Russell et al. (1996) estimate the 
rate of increase of stratospheric fluorine to be 6.6 per cent per year. 

The second possible reason for the scatter in Figure 22 may be traced 
to the shorter chemical lifetime of CH 4 in older air containing smaller amounts of 
CH 4 (less than 1 ppm). Figure 23 shows the chemical lifetimes of CH 4 and HF in 
tropical air calculated by the Max Planck model. The height interval in the plot (30 
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to 55 km) was chosen since it is here that the mixing ratios of HF and CH 4 most 
closely match those observed at lower altitudes close to the poles. It is thus probable 
that the vortex air we focus on in our study was once carried aloft to this region 
of rapid chemistry before being tranported poleward in the typical Brewer-Dobson 
circulation. The Figure shows the chemical lifetime of CH 4 declining rapidly with 
height above about 35 km. At these high altitudes over the tropics, meridional 
mixing may compete with chemical processes, and the CH 4 content of a parcel may 
become more sensitive to the path the parcel takes in the atmosphere. This increased 
sensitivity may lead to parcels of similar age close to the pole containing variable 
amounts of CH 4 . 

The third possible reason for the scatter in Figure 22 may have to do with 
the small-scale mixing processes within the polar night jet a s described by Pierce 
et al. (1994b). Luo et al. (1995) have noted that the slopes of scatterplots of 
CH 4 versus HF within the vortex differ from those of similar plots of these tracers 
outside the vortex. Small-scale mixing could bring together mid-latitude and vortex 
air and blur the tight linear correlation between these two tracers observed at other 
latitudes. 

Because of the scatter in Figure 22, we have chosen to sort parcels by 

c 

both CH 4 and HF content in order to accomplish our goal of identifying a series of 
parcels with similar age and temperature history. In other words, we have set limits 
on both the CH 4 and the HF content for those parcels we wish to examine most 
closely. For CH 4 we have chosen the interval of 0.75 to 0.9 ppm; for HF, 560 to 650 
ppt. In this way, the parcels represented by the outliers in Figure 22 are excluded 
from the analysis. The box in Figure 22 encloses those parcels we have chosen for 
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our study. These parcels contain mixing ratios of HF and CH 4 that are typical 
of air just poleward of the vortex edge (Pierce et al., 1994; Luo et ah, 1995). The 
chosen range of HF mixing ratios, well below the 900 ppt limit, excludes air that has 
experienced the rapid rates of NO photolysis mentioned above. We must emphasize 
that parcels characterized by CH 4 and HF mixing ratios outside the chosen ranges 
have also shown evidence of PSC-processing and ozone loss (e.g., Tuck et ah, 1995). 
These parcels were omitted from our study simply because during the time period 
under investigation HALOE did not observe a sufficient number of such parcels to 
unambiguously define trends in the evolution of short-lived trace gases. 

It is important to check whether parcels selected by this method show any 
trend in tracer mixing ratio over time, even within the small intervals chosen. Such 
a trend could bias the results discussed in the next section. Figure 24 shows a plot of 
HF versus day of year for October 1-23, 1992, on the 480 K surface for those parcels 
specified by our constraints on HF and CH 4 mixing ratios. The dotted line shows 
the average latitude of the parcels selected on each day. The correlation coefficient 
R of mixing ratios with day is only 0.02. Figure 25 shows a similar plot for CH 4 
using the same constraints. Again no trend over time is observed ( R = -0.04). 

6.3. Observations of Chemical Change 

To study the chemistry of air parcels close to the pole during Antarctic 
spring, HALOE observations during October of 1992 were interpolated to surfaces of 
constant potential temperature. Following the reasoning of the preceding section, we 
selected those parcels on the 480-K surface whose CH 4 mixing ratios fell between 0.75 



97 


and 0.9 ppm and whose HF mixing ratios fell between 560 and 650 ppt. Figure 26 
shows the ozone data for those parcels during the same time period as Figures 24 
and 25. Most of the points cluster below about 1 ppm although there are outliers. 

A closer look at the parcels with higher ozone shows that many of them 
are located in regions of steep horizontal gradients in HF or CH 4 mixing ratios. 
Figure 27 shows the distribution of HF on the 480-K surface during October 14-23, 
1992. During this time period, HALOE sunset observations gradually circled away 
from the vortex region. The triangles represent the parcels whose HF and CH 4 
mixing ratios are within the specified ranges. The empty triangles denote those 
parcels with O3 less than 1 ppm; the solid triangles, those parcels with O3 more 
than 1 ppm. Nine out of 10 parcels with high O3 are located at about 77 S between 
90 E and 170 W, a region where the HF distribution changes rapidly with latitude. 
Bithell et al. (1994) have noted that the HALOE observations in this region in 
October 1992 were made close to the strongest winds (i.e., within the jet core), 
where mixing between mid-latitude and vortex air was most likely to occur. Thus 
the high O3 mixing ratios are likely to be evidence of such mixing. Despite this 
mixing, the median of all the O3 points in Figure 26 is only 0.6 ppm. 

Figure 28 shows HC1 during the same time period. For this and subsequent 
data plots, we show the daily average of observations for those parcels that meet our 
criteria for CH 4 and HF and also contain ozone mixing ratios of less than 1 ppm. 
The plot shows a steep rise in HC1, from about 0.8 ppb to 2.2 ppb, during the first 
ten days of the period, followed by a leveling off to about 2.3 ppb. This trend is 
consistent with the increase in HC1 column amounts observed by Liu et al. (1992) 
in late October at McMurdo Station. 
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Some care must be taken in presenting the NO data since the sequence 
of measurements begin with sunrise observations and end with sunset. Because of 
nighttime odd nitrogen chemistry, sunset measurements of NO are not expected to 
match sunrise measurements the following morning, and comparing sunrise measure- 
ments from early October to the sunset measurements from the rest of the month 
may be misleading. Immediately after sunset, when N0 2 photolysis ceases, all NO 
is converted rapidly to N0 2 . During the night N0 2 converts more slowly to N 2 Os. 
In regions of high CIO, the nighttime chemistry must also include the three-body 
reaction of CIO + N0 2 to make C10N0 2 . At dawn, N0 2 photolysis begins renewed 
production of NO, but the sunrise sum of N0-fN0 2 , or NOx, will be less than 
the sunset sum because of the nighttime conversion N0 2 to N 2 0s and, possibly, to 

C10N0 2 . 

To correct the sunset values of NO to sunrise values, we first assume that 
CIO abundances are small by the time HALOE begins making sunset measurements, 
on October 9. As is discussed below, it is doubtful that much CIO remains by this 
point. We then follow the procedure outlined in Section 2.4.2 to calculate sunrise 
N0+N0 2 from the N0+N0 2 of the previous sunset. In Figure 29, the daily average 
NO for sunrise (October 1 through 7) and corrected sunrise (October 10 through 
23) is plotted for the same set of parcels as was selected for Figure 28. Low levels 
of NO are observed for the first few days, followed by a rapid rise in mixing ratio. 
Correcting the sunset NO to sunrise value decreases the mixing ratio by an average 
of 15 per cent. Even with the correction, the NO mixing ratios increase by about a 
factor five for the time period. 
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Figure 30 shows the daily average of HALOE measurements of N0 2 for the 
same set of parcels. As for NO, the N0 2 data in the second part of the period were 
corrected to show sunrise mixing ratio on the following morning. Again, making 
the correction decreases the N0 2 values by about 15 per cent. The trend in N0 2 
resembles that of NO with, however, much more scatter, perhaps due to the greater 
difficulty in measuring N0 2 at low altitudes (Gordley et al., 1996). Nonetheless a 
general increase in mixing ratio is apparent from values of about 0.15 ppb to values 
of more than 0.4 ppb. 


6.4. Model Results 

To interpret our observations of trace species during Antarctic spring, cal- 
culations were performed with the photochemical box model descibed in Chapter 3. 
The model was initialized at dawn in the following manner. First, since little work 
has been done to infer inorganic chlorine from either CH 4 or HF abundances, we 
used the method of Woodbridge et al. (1995), which infers Cly from N 2 0. In this 
method, total chlorine — the sum of both inorganic and organic chlorine — is first 
calculated using N 2 0 observations, the 1992 tropospheric mixing ratio of total chlo- 
rine (3670 ppt), and a growth rate of tropospheric chlorine of 2.8 per cent (WMO 
1992). An empirical relationship between organic chlorine and N 2 0 is then imple- 
mented to determine the inorganic fraction of total chlorine. In our analysis, we 
used N 2 0 observations made by the Cryogenic Limb Array Etalon Spectrometer 
(CLAES), another instrument on UARS, 13 days prior to the arrival of HALOE in 
the region (Roche et al., 1993). In the same regions that HALOE later measured 
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HF and CH 4 values between the limits specified above, N 2 0 observations averaged 
about 130 ppb and inferred Cly mixing ratio, about 2.7 ppb. We chose 2.7 ppb as 
our input value for Cly. As an aside, the value of N 2 0 mixing ratio obtained from 
CLAES data for the specified range of CH 4 agrees well with observed correlations 
of CH 4 and N 2 0 from ER -2 data (e.g., Kelly et al., 1989). 

We partitioned Cly into three species: HC1, C10N0 2 , and CIO. Based on 
HALOE data, a value of 1.0 ppb was assigned to HC1. Although C10N0 2 was also 
measured by CLAES, the lifetime of CIONO 2 at this time of year is on the order of 
days, and chlorine partitioning may well have changed in the days before HALOE 
arrived in the region. To determine an input value for C10N0 2 , we first note that 
the sum of HALOE N0+N0 2 increased about 1.0 ppb during the period. Our 
calculations show that little HNO 3 can photolyze to make NO+NO 2 during the 23 
days under study. We have therefore chosen to assign the missing nitrogen, 1.0 ppb, 
to CIONO 2 and to make 1.0 ppb the input value of that gas. (Our calculations show 
that nearly all the C10N0 2 photolyzes during the time period, so it is doubtful that 
the initial mixing ratio exceeded 1.0 ppb.) That choice leaves 0.7 ppb of remaining 
Cly, which we have assigned entirely to CIO. Given these initial values for CIO and 
CIONO 2 , the model calculates the mixing ratios of Cl and CI 2 O 2 at each timestep. 

For the NOy species, we again relied on HALOE data for the input values of 
NO and N0 2 ( 0.1 ppb each). For HNO 3 , we looked at CLAES data from the previous 
month in the same region that HF and CH 4 mixing ratios met our specifications, and 
found that HNO 3 mixing ratios there averaged about 3.2 ppb. Total active nitrogen 
(CIONO 2 + HNO 3 + NO + NO 2 ) was thus input as 4.4 ppb, a value consistent 
with previous observations of denitrification in the region a few degrees poleward of 
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the vortex edge (Fahey et al., 1990). Input mixing ratios of 0 3 and CH 4 were both 
taken from HALOE measurements — 0.6 ppm in the case of 0 3 and 0.8 for CH 4 . 
Aerosol surface areas were input as 3.6 /z 2 /cm 3 , or half the value of that calculated 
from HALOE aerosol extinction data (M. Hervig, personal communication). Only 
half the calculated value was used since comparisons with balloon data show some 
discrepancies, with HALOE overestimating aerosol surfaces areas (Hervig et al., 
1996). 

The model was run for 23 days at 20 km and 72 S, with increasing sunlight 
on each day. Temperature was held fixed at 206 K, which was the average tem- 
perature for the time interval for those parcels with less than 1 ppm of ozone (i.e., 
those parcels that had experienced the least mixing with mid-latitude air). The 0 3 
column amount, which influences photolysis rates, was also held fixed at 3.8 xlO 18 
molecules/cm 2 , the average HALOE value for the period. Model results show that 
0 3 declined slightly (about 0.5 ppm) during the first ten days, but then held con- 
stant at about 0.55 ppm for the remaining days. Figure 31 shows the evolution of 
C10N0 2 , CIO, and HC1 during the 23 days of the model run. While C10N0 2 val- 
ues are high in the first part of the period, they fall rapidly as the HC1 abundance 
quickly builds. By the 15th d^r, HC1 has reached a value of 2.6 ppb, somewhat 
higher than the 2.3 ppb observed by HALOE. This represents a model overestimate 
of about 13 per cent. 

Figure 32 presents model results for four active nitrogen species — HN0 3 , 
NO 2 , NO, and CIONO 2 . The trace gas C10N0 2 is shown a second time for com- 
parison with other nitrogen-containing species, and NO and N0 2 are shown here 
as the sum N0-t-N0 2 . Calculated HN0 3 shows an initial decline during the period 
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due to increased photolysis as the days lengthen. By about the 15th day, however, 
when much of the chlorine in C10N0 2 has shifted to HC1, the nitrogen contained in 
CIONO 2 has shifted mainly to N0+N0 2 . Some of the nitrogen ends up in HNO3, 
which returns to its initial value. The abundance of N0-f-N0 2 rises steadily over 
the course the period, reaching a level of about 1.0 ppb, a value about the same as 
the observations. 


6.5. Discussion 

Using HF and CH 4 as tracers to identify parcels of similar air, we have 
detected a rapid rise of HC1 during October in air that has experienced springtime 
vortex conditions in the lower stratosphere over Antarctica. The rise in HC1 is ac- 
companied by similar rises in NO and N0 2 . Ozone mixing ratios show no trend 
during the time interval. Following the reasoning laid out in Section 3, we can as- 
sume that Cly and NOy levels remain roughly constant for all the parcels considered 
in this study. Thus, the observed changes in the short-lived species are probably not 
due to mixing processes bringing in air of higher Cly or NOy content to the region, 
but to rapid chemistry occurring in the aftermath of extensive ozone loss. The model 
results support our conclusion tjjat these changes are indeed due to chemistry. 

Douglass et al. (1995), using MPV as a tracer on the 460 K surface, saw 
similar rises in the three trace gases under investigation here — HC1, NO, and 
N0 2 . While the trends detected by both methods show some scatter, those using 
tracer analysis show less, especially for HC1 and NO. In addition, our model results 
better matched observations than the predictions of Douglass et al. (1995), probably 
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because our input values for Cly and NOy were signficantly lower. In the Douglass 
et al. study, model results overpredicted these gases in the latter part of October by 
about 25 per cent for HC1 and more than 50 per cent for the sum of NO+NO2. 

To understand the dramatic chemical changes occurring in the partitioning 
of the Cly family, one must look closely at the set of reactions that govern the 
partitioning of Cly. HC1 is formed by the following reaction 


Cl + CH 4 -► HC1 + CH 3 


( 6 . 1 ) 


and destroyed by 


HC1 + OH Cl +H 2 0 (6.2) 

Both reactions occur slowly in the lower stratosphere, and the lifetime of HC1 in 
this region is on the order of weeks. Even with abundant CH 4 and low levels of 
O3, Reaction 6.1 takes place several orders of magnitude more slowly than the 
corresponding reaction of atomic chlorine with 0 3 : 


Cl + 0 3 — * CIO *f- 0 2 (6.3) 

In order for HC1 abundances to rise rapidly in springtime Antarctica, the Cl 
fraction of available chlorine mu5t increase at the expense of other reactive chlorine 
species. As O3 levels fall, the rate of Reaction 6.3 declines slightly, thus increasing 
Cl abundances and the rate of Reaction 6.1. An even more effective mechanism to 
increase the Cl fraction of Cly involves two competing reactions of the NO radical: 

NO +C10 Cl +N0 2 


NO T 0 3 


O2 T NO2 


(6.4) 

(6.5) 
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Under normal stratospheric conditions of abundant ozone and little CIO, Reaction 
6.4 occurs many times slower than Reaction 6.5. Under the perturbed conditions 
that are typical of springtime over Antarctica, with low 0 3 and high CIO, the rate 
of Reaction 6.4 may become similar to, or faster than, the rate of Reaction 6.5. 

Figure 33 presents the rates of three reactions important to springtime 
chemistry in the stratosphere over Antarctica. The plot shows the daily average rate 
of Reactions 6.1, 6.4, and 6.5 for each day of the period of the HALOE observations. 
As in the preceding model plots, the altitude of the calculations is 20 km and the 
latitude is 72 S. The model has been initialized in the same manner as described 
above. For clarity, the rate of Reaction 6.1 has been multiplied by a factor of .01, 
while the other two reaction rates have been multiplied by 10~ 5 . In the early days 
shown in the plot, the rates of all the reactions considered are increasing with time. 
After about the 9th day, however, as the rate of the Reaction 6.5 begins to overtake 
the rate of the competing Reaction 6.4, the rate of HC1 production peaks and then 
declines. 

The plot suggests that it is the relative magnitudes of the rates of two NO 
reactions, Reactions 6.4 and 6.5, that determine the pace of HC1 recovery during 
Antarctic spring. In the early days of the simulation, when large abundances of CIO 
are present, Reaction 6.4 proceeds more rapidly than Reactions 6.5, and chlorine 
as a consequence is shifted from CIO to HC1 via atomic chlorine. With less CIO in 
abundance, CIONO 2 concentrations also fall, due to a decreased rate of the following 
reaction: 


CIO +N0 2 +M 


CIONO2 +M. 


( 6 . 6 ) 
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Less C10N0 2 , in turn, means that more odd nitrogen is available in the form of NO 
and NOj, and Reaction 6.4 may proceed even more quickly in a positive feedback 
loop. At some point, however, so much chlorine has been shifted to HC1 that 
Reaction 6.4 slows down. The rate of Reaction 6.5, on the other hand, continues to 
increase steadily as NO abundances increase, and will eventually overtake the rate 
of Reaction 6.4. 

To test these ideas, we can run the box model with varying amounts of 
initial ozone abundance. Figure 34 shows the evolution of chlorine species for an 
initial O3 of only 0.3 ppm, or half the initial O3 of the calculations described above. 
Figure 35 shows the same set of reactions as Figure 33, but with the smaller ozone 
input value. With less initial ozone, HC1 production peaks earlier, on about day 5, 
but the magnitude of the peak is about double that with 0.6 ppm of ozone. Again, 
the rate of HC1 production declines once the rate of Reaction 6.5 overtakes that of 
Reaction 6.4. The earlier peak of HC1 production is reflected in Figure 34, with a 
surge of HC1 produced at the expense of C10N0 2 early in the time period. 

In Figure 36, we show the model calculation of Cly with an initial input of 
2.4 ppm of ozone, or 4 times as much ozone as in the first model run. Figure 37 shows 
the same set of relevant reactions as above. With 2.4 ppm of ozone, HC1 mixing 
ratios climb much more slowly, gaining only about 0.1 ppb during the entire period. 
As the greater ozone abundance shifts Cl to CIO and NO to N0 2 , C10N0 2 mixing 
ratios climb to a daytime value of about 1.5 ppb. CIO levels decline somewhat but 
still reach noontime values greater than 0.1 ppb throughout the period. In Figure 
37, the rates of all three key reactions rise steadily during the 23 days. However, 
the rates of Reactions 6.1 and 6.4 are much slower than the corresponding rates in 
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the previous model runs, and even by the end of the period have not yet reached 
the same values as the peak values obtained before. Reaction 6.5 rises most rapidly 
of the three reactions, and by day 15 has overtaken Reaction 6.4. 

As Douglass et al. (1995) and others have shown, one may express the 
balance of Cl to CIO mathematically by assuming that all of the relevant radicals are 
in steady state, and that the reactions of greatest importance under these conditions 
are Reactions 6.3, 6.4, 6.5, and the photolysis of N0 2 : 


N0 2 +hu -* NO +0. 


(6.7) 


With these assumptions, the ratios of N0/N0 2 and Cl/CIO may be written: 

[NO ] _ ^ 6.7 

[N0 2 ] “ MOs] + jmcio] 


[Cl] = ^ 6.4 [NO] 
[CIO] ^ 6 . 3 ( 03 ] 


(6.9) 


Solving for NO in Equation 6.8 and rewriting Equation 6.9, we have: 


[Cl] = fc 6 , 4 J 6 , 7 [ NQ 2 ] 

[CIO] h. 5 k 6 . 3 [0 3 ) 2 + k 6A k 6 . 3 {C\0][0 3 ] (6,10 ' 

Thus, the steady-state ratio of Cl to CIO depends on the inverse of the square of 

* 

the O3 abundance. In regions of low ozone, chlorine is shifted rapidly from CIO to 
Cl, and then on to HC1. According to the model, nearly all the available chlorine 
(2.6 out of 2.7 ppb) has been shifted to HC1 by the end of the 23-day period, while 
observations show that only about 2.3 ppb of HC1 has been formed by the 23rd 
day. Reasons for this minor discrepancy are not clear. In any case, since HC1 is so 
long-lived, the repartitioning of the Cly family quickly shuts off further 0 3 loss due 
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to chlorine catalysis. As a final test of the effect of the reaction of CH 4 with atomic 
chlorine, we have calculated the evolution of O3 during this time period in a run 
that excludes Reaction 6.1 from the reaction scheme. The plot also shows model 
ozone for the standard run. As noted above, ozone in the standard run drops only 
slightly in the first few model days and then levels off. Model ozone in the run that 
neglects Reaction 6.1 shows a precipitous drop to about two-thirds its initial value. 

6.6. Summary 

HALOE observations of trace species within the Antarctic vortex in late 
September and through October of 1992 suggest a rapid repartitioning of strato- 
spheric chlorine from reactive species to the reservoir species HC1. This dramatic 
chemical change may be detected by using HF and CH 4 together as tracers to se- 
lect a series of air parcels of similar photochemical age, having similar abundances 
of Cly and NOy and similar temperature histories. The rapid repartitioning of 
chlorine comes about when ozone abundances are very low and abundances of re- 
active chlorine are very high. Under these conditions, (1) the rate of the reaction 
Cl "t~ O3 ♦ CIO -f O2 slows, and (2) the rate of the reaction NO + CIO — +■ NO2 + Cl 
becomes faster than the rate of the competing reaction NO + O3— ► NO2 + 0 2 . 

v 

Both consequences lead to increases in Cl at the expense of CIO. Increased abun- 
dances of Cl allow chlorine to be shifted much more rapidly to HC1 through the 
reaction Cl -f CH 4 — > HC1 + CH3. Since HCl has a long lifetime, on the order of 
weeks, the sudden rise in HCl effectively shuts down the chlorine catalytic processes 
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that deplete ozone. The shift in chlorine may have importance as filaments of vortex 
air peel off the vortex edge in spring and are transported toward mid-latitudes. 


* 
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Figure 21. Polar stereographic contour plot of HALOE sunrise HF on the 480- K 
surface from September 7 to October 5, 1992. Contours are in parts per billion. 
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Figure 22. Correlation of HALOE CH4 vs HF mixing ratios between 15 and 25 km 
and 62 and 78 S during the time period October 1-23, 1992. 
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Figure 23. Chemical lifetimes of HF and CH 4 , calculated by a two-dimensional, pho- 
tochemical model. The circles denote HF production lifetime due to CFC photolysis; 
the triangles show net CH 4 loss by reaction with 0( l D), Cl, and OH. 
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Figure 24. HALOE sunrise and sunset HF observations (stars) on the 480 K surface 
for October 1-23, 1992. Only those parcels containing 560 to 650 ppt of HF and 0.75 
to 0.9 ppm of CH 4 are shown. The dotted line shows the average daily latitudes 
observed for the parcels within the specified limits of HF and CH 4 . 
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Figure 25. HALOE sunrise and sunset CH 4 observations, as in Figure 24. 
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Sunrise and sunset 03 at 480 K 



Figure 26. HALOE sunrise and sunset O 3 observations on the 480-K surface. Only 
those parcels containing 560 to 650 ppt of HF and 0.75 to 0.9 ppm of CH 4 ppt are 
shown. 
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Sunset HF at 480 K 
from 6 Oct 1992 to 14 Oct 1992 
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Figure 27. Polar stereographic contour plot of HALOE sunset observations of HF 
on the 480 K surface from October 14 to October 31, 1992. The outer latitude ring 
is 70 S. The triangles represent the parcels with a HF and CH 4 content within the 
specified limits (see text). The solid triangles denote those parcels containing more 
than 1 ppm of ozone; the unfilled triangles, those parcels containing less than 1 ppm. 
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Sunrise and sunset HC1 at 480 K 



Figure 28. Daily average HALOE sunrise and sunset HC1 at 480 K for the same 
time period and HF and CH 4 constraints as Figure 26. 
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Sunrise and sunset NO at 480 K 



Figure 29. Daily average HALOE NO at 480 K for the same time period and HF 
and CH 4 constraints as in Figure 26. The data points before October 10 represent 
sunrise observations; after that date the data have been corrected to sunrise. 


Latitude 




118 


Sunrise and sunset N02 at 480 K 



Figure 30. Daily average HALOE NO 2 at 480 K for the same time period and HF 
and CH 4 constraints as in Figure 26. The data points before October 10 represent 
sunrise observations; after that date the data have been corrected to sunrise. 
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Figure 31. Calculated HC1, CIONO 2 , and CIO at 20 km during the same time period 
as the HALOE observations. 
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Figure 32. Calculated HNO3, CIONO2, and NO+NO2 at 20 km during the same 
time period as the HALOE observations. 



121 



Day in October 


Figure 33. Rates of three key reactions in springtime chemistry over Antarctica 
calculated with the same model conditions as above. The rate of the Reaction Cl 
+ CH4 has been multiplied by a a factor of .01, while the other two reaction rates 
have been multiplied by 10 -5 . 



Cly Species (ppb) 


122 



1 5 <*9 13 17 21 

Day in October 


Figure 34. Calculated HC1, CIONO 2 , and CIO at 20 km as Figure 31, but with only 
0.3 ppm initial ozone. 
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Figure 35. Rates of key reactions in springtime chemistry as in Figure 33, but with 
only 0.3 ppm initial ozone. 
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Day in October 


Figure 36. Calculated HC1, CIONO 2 , and CIO at 20 km as in Figure 31, but with 
2.4 ppm initial ozone. 
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Figure 37. Rates of key reactions in springtime chemistry as in Figure 33, but with 
2.4 ppm initial ozone. 





CHAPTER 7 


CONCLUSIONS. 

This work shows the utility of tracer analysis as a means to distinguish 
between change due to transport processes in the stratosphere and change due to 
chemistry. Tracer material serves as a kind of clock that assigns a photochemical age 
to the air parcels observed. We may assume that parcels with similar abundances of 
tracer material have been exposed to similar quantities of solar radiation and thus 
contain similar abundances of the chemical families that derive from source gases, 
such as N 2 O or the CFCs. Thus for any series of air parcels, observations of just one 
or two member species of a family provide sufficient information to draw conclusions 
about the partitioning of that family over timescales of weeks or even years. 

In two studies, we used tracer analysis to follow chemical processes of 
importance to the stratospheric ozone budget. In the first study, we examined the 
response of NOx and ozone at 17 mb to the presence of aerosol leftover from the 
1991 eruption of Mount Pinatubo. One important conclusion of this study was that 
NOx in this region of the stratosphere plays a dominant role in ozone loss chemistry. 
This result could have implications in the event that a fleet of high-speed civilian 
transport aircraft is launched, and aircraft emissions of reactive nitrogen reach this 
high in the stratosphere. In the second study, we tracked the abundances of chlorine 
and nitrogen species in the wake of significant ozone depletion over Antartica in 
spring. Here we found that under these conditions, a dramatic shift in chlorine from 
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radical species to the reservoir species HC1 takes place, shutting down further ozone 
loss via catalytic chlorine cycles. 
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